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Abstract 

Clear-air  turbulence  (CAT)  prediction  is  vitally  important  to  military  aviation  and 
the  successful  completion  of  Department  of  Defense  (DoD)  operations  such  as  air  to  air 
refueling  and  new  national  defensive  weapon  systems  such  as  directed  energy  platforms. 
The  unique  mission  requirements  of  military  aircraft  often  require  strict  avoidance  of 
turbulent  regions.  Traditionally,  weather  forecasters  have  found  it  difficult  to  accurately 
predict  CAT.  In  order  to  forecast  regions  where  CAT  might  occur,  forecasters  must  first 
determine  the  location  of  breaking  waves  caused  by  either  Kelvin-Helmholtz  instabilities 
or  topographically  forced  internal  gravity  waves  (mountain  waves)  in  the  atmosphere. 

The  United  States  Air  Force  (USAF)  15*  Operational  Weather  Squadron  (15th  OWS) 
requested  an  updated  method  of  predicting  CAT  and  this  request  was  ranked  as  one  of  the 
highest  priority  research  needs  by  the  HQ  USAF  Director  of  Weather,  Deputy  Chief  of 
Staff  for  Air  and  Space  Operations. 

A  new  method  of  forecasting  turbulence  was  developed  in  this  work  and  the 
operational  model  was  delivered  to  the  15*  OWS  for  immediate  inclusion  into  their 
operations.  This  method  combines  output  from  the  Knapp-Fllrod  index  and  the  Naval 
Research  Faboratory’s  Mountain  Wave  Forecast  Model  (MWFM)  onto  a  single  chart. 
Displaying  these  tools  together  allows  forecasters  to  view  both  causes  of  CAT 
simultaneously.  Furthermore,  a  new  visualization  tool  is  developed  that  allows  a 
forecaster  to  view  several  layers  at  the  same  time  as  well  as  a  composite  chart  to  greatly 
reduce  the  time  required  to  produce  turbulence  charts  by  OWS  forecasting  centers 


IV 


worldwide.  Tests  of  forecast  accuracy,  as  determined  by  pilot  reports  (PIREPS),  between 
charts  currently  produced  by  USAF  OWSs  and  this  new  method  were  compared,  with  the 
new  method  producing  far  superior  forecast  results.  This  method  revolutionizes  the  way 
USAF  OWSs  produce  turbulence  charts  to  help  the  Air  Force  satisfy  each  of  its  core 
competencies;  especially  information  superiority,  rapid  global  mobility,  and  precision 
engagement. 
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AN  AUTOMATED  METHOD  OF  PREDICTING  CLEAR-AIR  TURBULENCE 


I.  Introduction 

Over  the  history  of  aviation,  atmospheric  turbulence  has  had  a  huge  impact  on 
flight  cost  and  safety.  For  example,  turbulence  cost  United  States  commercial  and 
military  aircraft  over  $25  million  annually  from  1963  to  1965  (Hopkins  1977).  More 
recently,  turbulence  related  incidents  accounted  for  over  half  of  all  the  commercial  carrier 
accident  and  incident  reports  between  the  years  1990  and  2000,  according  to  the  National 
Transportation  Safety  Board  (Sharman  et  al.  2000).  Generally,  there  are  two  categories 
of  atmospheric  turbulence  of  concern  to  aviators:  turbulence  associated  with  convective 
clouds  and  turbulence  that  occurs  outside  of  convective  clouds.  Pilots  and  weather 
forecasters  are  primarily  concerned  with  turbulence  occurring  outside  of  convective 
clouds,  or  clear-air  turbulence  (CAT),  as  it  rarely  offers  any  visual  indications. 

CAT  is  defined  as  any  turbulence  in  the  free  atmosphere  of  interest  to  aerospace 
operations  that  is  not  in  or  adjacent  to  visible  convective  activity,  including  turbulence 
found  in  cirrus  clouds  not  in  or  adjacent  to  visible  convective  activity  (Hopkins  1977). 
There  are  two  widely  accepted  causes  of  CAT;  either  standing  waves  in  the  lee  of  a 
mountain  barrier  (mountain  waves)  or  vertical  wind  shear  in  a  stably  stratified  layer 
(Kelvin-Helmholtz  instabilities)  (Hopkins  1977).  Current  methods  of  prediction  have  a 
difficult  time  accurately  forecasting  CAT.  This  forecast  challenge  exists  in  part  because 
forecasters  must  use  synoptic  scale  analysis  techniques  and  coarse  resolution  numerical 
models  to  predict  a  mesoscale  phenomenon  (Ellrod  1989). 


1 


1. 1  Statement  of  Problem 


Whether  supporting  commereial  or  military  aircraft,  forecasters  spend  significant 
time  and  effort  predicting  CAT.  For  example,  USAF  weather  personnel  make  turbulence 
charts  for  12,  18,  24,  30,  36,  and  48  hour  forecasts  once  per  day  to  support  United  States 
Army  and  Air  Force  aviation.  These  charts  should  include  regions  of  light,  moderate, 
severe,  and  extreme  turbulence  time  phased  for  three  hours  before  and  after  the  forecast 
hour  between  10,000  and  50,000  feet  mean  sea  level  (MSL)  (Air  Force  Manual  15-129 
2001).  However,  in  practice,  most  weather  squadrons  only  forecast  for  moderate  or 
greater  (MOG)  turbulence.  The  range  between  10,000  and  50,000  feet  is  also  represented 
as  FL100-FL500.  Past  experience  has  shown  these  charts  take  significant  time  to 
prepare. 

USAF  weather  forecasters  use  numerous  methods,  some  of  which  are  subjective, 
to  predict  CAT  (15  OWS  2003).  Subjective  methods,  such  as  conceptual  models  and 
rules  of  thumb,  rely  heavily  on  individual  forecaster  skill  and  experience,  both  of  which 
vary  greatly  between  USAF  weather  forecasters.  Since  USAF  weather  centers  predict 
turbulence  for  large,  geographically  diverse  portions  of  the  earth,  subjective  methods  are 
largely  responsible  for  the  substantial  length  of  time  required  in  building  these  charts. 

Further  increasing  the  time  required  to  build  turbulence  charts,  different  methods 
of  CAT  prediction  capture  one  cause  of  turbulence  or  the  other,  but  not  both.  This 
limitation  forces  forecasters  to  use  multiple  tools  while  building  their  products.  The 
method  developed  here  provides  a  tool  that  displays  indices  representing  both  causes  of 
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CAT  on  the  same  graphic,  allowing  the  forecaster  to  quickly  determine  regions  where 
either  might  occur. 

1.2  Research  Objectives 

The  primary  goal  of  this  research  is  to  provide  Air  Force  weather  squadrons  with 
an  invaluable  tool  for  timely  and  accurate  prediction  of  CAT.  Specifically,  the  objectives 
of  this  research  are  to  provide  USAF  weather  personnel  with: 

1 .  A  method  to  quickly  determine  likely  regions  of  turbulence  caused  by  either 
Kelvin-Helmholtz  instabilities  or  mountain  waves; 

2.  An  objective  tool  requiring  little  human  interpretation  that  provides  superior 
forecasts  to  those  currently  requiring  significant  human  interpretation; 

3.  A  web-based  visualization  tool  that  provides  a  three-dimensional,  layer-by-layer 
view  of  the  atmosphere;  and, 

4.  An  initial  determination  of  the  potential  benefit  of  using  the  Mountain  Wave 
Forecast  Model  (MWFM)  to  predict  turbulence  in  the  troposphere. 

1.2  Research  Approach 

The  Ellrod-2  index,  as  described  by  Ellrod  and  Knapp  (1991),  was  chosen  to 
depict  Kelvin-Helmholtz  instabilities.  It  computes  a  value  based  on  vertical  wind  shear, 
horizontal  deformation,  and  convergence.  This  index  is  built  within  the  Operational 
Production  System  Version  2  (OPS  II)  software  used  by  most  OWSs.  Eurther 
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substantiating  its  use,  this  index  eompared  favorably  with  other  indiees  in  a  study  by 
Brown  et  al.  (2000). 

The  MWFM  version  1  and  2  were  ehosen  to  prediet  mountain  waves  in  this  study. 
These  models  were  developed  at  the  Naval  Researeh  Laboratory  (NRL)  to  prediet  regions 
of  stratospherie  turbulenee  eaused  by  breaking  mountain  waves.  The  MWFM  is 
eomputed  from  atmospheric  modeled  data  and  a  global  ridge  database.  Application  of 
the  model  to  two  case  studies  of  tropospheric  turbulence  suggested  it  might  apply  there  as 
well  (Eckermann  et  al.  2000).  In  this  study,  both  versions  of  the  model  were  tested  to 
determine  if  both  indices  added  value  to  a  forecast. 

The  indices  were  calculated  over  the  continental  United  States  (CONUS)  using 
the  National  Center  for  Environmental  Prediction’s  (NCEP)  Global  Eorecast  System 
(GPS)  model,  formerly  known  as  the  Aviation  model  (AVN).  The  MWPM  retrieved  the 
data  using  file  transfer  protocol  (PTP)  and  stored  the  data  locally  for  access  by  a  program 
that  computes  the  Ellrod-2  index.  To  facilitate  timely  completion  of  this  research,  only 
the  analysis,  12-hour  forecast,  and  24-hour  forecast  were  produced. 

Initially,  the  output  from  the  indices  were  compared  with  pilot  reports  (PIREPS), 
as  provided  by  the  Air  Porce  Combat  Climatology  Center  (AECCC),  to  determine  the 
threshold  for  each  index.  After  determining  the  thresholds,  the  indices  were  combined 
into  a  single  graphic  for  each  50  mb  layer  of  the  atmosphere  that  contained  PE100-PE500 
and  integrated  into  a  multiple  layer  viewer  for  forecasting.  Additionally,  a  chart  was 
made  containing  a  composite  view  of  all  layers.  Pinally,  14  days  of  forecasts  using  these 
charts  were  developed.  The  accuracy  of  these  forecasts  and  forecasts  produced  by  the 
15*  Operational  Weather  Squadron  (OWS)  were  verified  against  PIREPS  data. 
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The  remainder  of  this  thesis  provides  details  on  the  indiees  ehosen  and  the 
methods  used  to  forecast  turbulence.  Chapter  two  provides  background  information  on 
the  causes  of  CAT  and  the  chosen  indices.  Then,  chapter  three  describes  the 
implementation  of  the  forecast  method.  Chapter  four  discusses  the  results  of  the 
threshold  determination  and  the  forecast  comparison.  Finally,  chapter  five  provides 
conclusions  and  suggests  future  research  regarding  this  topic. 
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II.  Literature  Review 


2.1  Kelvin-Helmholtz  Instability 

Meteorologists  have  long  reeognized  and  aeeepted  breaking  waves  eaused  by 
Kelvin  Helmholtz  instabilities  as  a  primary  eause  of  elear-air  turbulenee.  Kelvin- 
Helmholtz  instabilities  oeeur  in  regions  of  strong  vertieal  wind  shear  (VWS)  in  stably 
stratified  air  (Hopkins  1977).  Vertical  wind  shear  is  simply  the  change  in  wind  speed 
with  height  while  stable  stratification  requires  more  explanation. 

Stably  stratified  air  suggests  relatively  small  vertical  density  differences.  A 
modified  form  of  the  ideal  gas  law  relates  pressure,  temperature,  and  density  to  each 
other  by: 

P  =  pRT  (1) 

where  P  is  the  pressure,  p  is  density,  R  is  the  dry  air  gas  constant  287.04  J  ■  Kg~^  ■  K”' , 
and  T  is  the  temperature.  Typically,  atmospheric  pressure  and  temperature  decrease  with 
height.  Equation  1  suggests  that  changing  pressure  or  temperature  while  holding  the 
other  constant  changes  the  density.  The  equation  suggests  decreasing  pressure  causes 
density  to  decrease  while  decreasing  temperature  causes  density  to  increase.  Because  the 
effects  from  decreasing  pressure  outweigh  the  effects  from  decreasing  temperature 
(pressure  decreases  more  rapidly  with  height),  the  overall  effect  is  that  density  decreases 
with  height.  However,  the  more  rapidly  the  temperature  decreases  for  a  given  pressure 
decrease,  the  more  gradual  the  decrease  in  density.  Essentially,  the  more  rapidly  the 
temperature  decreases  with  height  the  greater  the  stable  stratification  of  the  air. 


To  simplify  the  mathematics  involved  in  determining  the  stability,  potential 
temperature  is  used.  Potential  temperature  is  defined  as  the  temperature  a  parcel  of  air 
would  have  if  it  were  compressed  or  expanded  adiabatically  from  its  existing  pressure 
and  temperature  to  a  reference  pressure  level  (Wallace  and  Hobbs  1977).  Basing 
calculations  on  potential  temperature  instead  of  standard  temperature  simplifies  the 
determination  of  the  stratification  of  the  atmosphere. 

Mathematically,  the  Richardson  number  (Ri)  relates  the  two  necessary  conditions 
for  Kelvin-Helmholtz  instabilities  as  a  ratio: 


Ri 


0  dz 

dv' 


dz 


(2) 


do 


where  g  is  the  acceleration  due  to  gravity,  0  is  the  potential  temperature,  —  is  the 

dz 

change  in  potential  temperature  with  height  (representing  the  stability  of  the  layer),  and 

dv 

- is  the  vertical  wind  shear.  In  order  for  turbulence  to  occur,  the  Richardson  number 

dz 

must  be  less  than  a  critical  value  of  0.25  (Thorpe  1969),  which  requires  some 
combination  of  low  static  stability  and  high  vertical  wind  shear.  It  is  hypothesized  that  a 
minimum  value  of  wind  shear  is  required  for  these  instabilities  to  form  regardless  of  the 
value  of  the  Richardson  number.  This  implies  that  meeting  the  0.25  threshold  is  a 
necessary,  not  a  sufficient,  condition  for  CAT  not  associated  with  mountain  waves 
(Ellrod  and  Knapp  1991). 
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It  is  also  important  to  note  that  Kelvin-Helmholtz  instabilities  are  possible  for 
small  wavelengths,  regardless  of  the  wind  shear  and  statie  stability  (Ray  1986). 

However,  small  wavelengths  do  not  eause  suffieient  turbulenee  to  influenee  aireraft. 
Aecording  to  Sharman  et  al.  (2000),  these  waves  need  to  be  on  the  scale  of  100  m  long  to 
affect  large  airframes. 

Often  these  waves  occur  in  air  with  insufficient  moisture  for  clouds  to  form. 

When  this  happens  they  do  not  have  any  visible  billow  clouds  associated  with  them.  If 
sufficient  moisture  is  present  these  waves  look  similar  to  breaking  ocean  waves.  Like  an 
ocean  wave,  the  larger  the  amplitude  of  the  wave  prior  to  it  breaking,  the  more  intense  the 
resulting  turbulence. 

Forecasters  have  relied  on  synoptic  pattern  recognition  and,  more  recently, 
numerical  model  indices  to  predict  these  instabilities,  with  limited  success.  Both  of  these 
approaches  attempt  to  use  large  scale  features  to  determine  regions  that  are  likely  to 
contain  these  relatively  small  scale  events  (Ellrod  1989).  Based  on  the  100  m  wavelength 
requirement,  numerical  model  grid  points  would  need  to  be  spaced  10  m  apart  to 
accurately  predict  these  waves  (Sharman  et  al.  2000).  Further  limiting  the  effectiveness 
of  these  indices  is  the  poor  temporal  resolution  of  both  synoptic  and  model  data.  Kelvin- 
Helmholtz  waves  are  short-lived  (often  less  than  30  minutes)  while  synoptic  and  model 
data  are  typically  valid  every  six  or  twelve  hours. 

With  the  introduction  of  computer  models,  researchers  have  recognized  the 
usefulness  of  indicators  other  than  the  Richardson  number  to  predict  CAT  associated  with 
Kelvin-Helmholtz  instabilities.  Although  a  Richardson  number  less  than  0.25  is 
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theoretically  a  good  indicator  of  turbulent  conditions,  it  is  not  well-suited  to  computer 
models  for  several  reasons: 

1 .  It  is  difficult  to  accurately  calculate  the  index  using  a  numerical  model  due 
to  coarse  vertical  resolution  (Ellrod  and  Knapp  1991). 

2.  Computations  using  data  from  rawinsondes,  the  best  source  of 
measurement,  often  lead  to  Richardson  values  between  0.25  and  1.0  in 
turbulent  regions  (Ellrod  and  Knapp  1991). 

3.  Values  of  Richardson  less  than  0.25  occur  in  both  turbulent  and  smooth 
areas  (Ellrod  and  Knapp  1991). 

4.  A  detailed  statistical  study  produced  poor  correlation  between  the 
Richardson  number  and  CAT  (Dutton  1980). 

Observed 


Eigure  1 .  Contingency  table  used  to  verify  turbulence  forecasts 

Recently,  Brown  et  al.  (2000)  conducted  a  comparative  study  of  1 1  CAT 
algorithms.  The  11  algorithms  tested  were  Brown-1  (Brown  1973),  Colson-Panofsky 
(Colson  and  Panofsky  1965),  DTP3,  DTP4,  DTP5  (Marroquin  1985,  1998),  Empirical 
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Dutton  (Dutton  1980),  Ellrod-1,  Ellrod-2  (Ellrod  and  Knapp  1991),  Endlich  (Endlich 
1964),  ITFA  (Sharman  et  al.  2000),  and  the  Richardson  number  (Drazin  and  Reid  1981; 
Kronebach  1964).  In  Brown  et  al.  (2000),  contingency  tables  such  as  the  one  depicted  in 
Figure  1  were  produced  for  each  index.  In  Figure  1,  observed  values  are  based  on 
PIREPS  while  forecasted  values  are  based  on  an  index. 


Table  1,  Statistics  used  in  index  studies 


Statistic 

Definition 

Description 

PODy 

A/(A+C) 

Probability  of  Detection  of  yes  observations. 
Probability  of  yes  observations  that  were 
correctly  forecasted. 

PODn 

D/(B+D) 

Probability  of  detection  of  no  observations. 
Probability  of  no  observations  that  were 
correctly  forecasted. 

FAR 

B/(A+B) 

False  Alarm  Rate.  Probability  of  yes 
forecasts  that  were  incorrect. 

CSl 

A/(A+B+C) 

Critical  success  index.  Number  of  correct 
yes  forecasts  relative  to  number  of  yes 
forecasts  or  observations. 

PC 

(A+D)/(A+B+C+D) 

Proportion  correct.  The  proportion  of  yes  or 
no  forecasts  that  were  forecasted  properly. 

Bias 

(A+B)/(A+C) 

Bias.  Frequency  of  yes  forecasts  relative  to 
yes  observations. 

HSS 

2(AD-BC)/ 

f(A+C)(C+D)+(A+B)(B+D) 

Heidke  Skill  Score.  PC  normalized  by  the 
number  of  forecasts  expected  to  be  correct. 

GSS 

(A-C2)/[(A-C2)+B+C) 
where  C2=(A+B)(A+C)/ 
(A+B+C+D) 

Gilbert  Skill  Score.  CSl  corrected  for  the 
value  of  YY  that  would  be  expected  by 
chance. 

TSS 

PODy+PODn-1 

True  Skill  Score.  Measure  of  discrimination. 

Using  PIREPS  and  automated  vertical  accelerometer  (AVAR)  observations  of  “g” 
forces,  as  reported  by  United  Airlines  flights,  as  observed  values  and  index  thresholds  as 
forecasted  values,  contingency  tables  were  produced  for  each  index  in  the  study.  The 
contingency  tables  were  used  to  compute  the  statistics  described  in  Table  1.  For  reasons 
discussed  in  detail  in  Section  2.3,  several  statistics  such  as  false  alarm  rate  (FAR),  bias. 
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and  critical  skill  index  (CSI)  can  be  somewhat  misleading.  When  considering  these 
statistics,  the  probability  of  detection,  yes  and  no  (PODy  and  PODn)  and  the  true  skill 
seore  (TSS)  are  fairly  reliable  indicators  of  an  index’s  forecast  ability. 


Table  2.  Comparison  of  indices  for  data  collected 


Algorithm 

Threshold 

PODy 

(All) 

PODy 

(MOG) 

POD„ 

(PIREPs) 

POD„ 

(AVARS) 

TSS 

AIRMETs 

~ 

0.57 

0.64 

0.70 

0.66 

0.34 

Brown-1 

0.58 

0.62 

0.63 

0.58 

0.25 

Colson-Panofsky 

-1000 

0.58 

0.63 

0.61 

0.53 

0.24 

DTF3 

0.7 

0.63 

0.68 

0.63 

0.63 

0.31 

DTF4 

2.5 

0.63 

0.67 

0.57 

0.60 

0.24 

DTF5 

0.15 

0.56 

0.59 

0.66 

0.68 

0.25 

Empirical  Dutton 

25 

0.62 

0.66 

0.57 

0.57 

0.23 

Ellrod-1 

0.61 

0.66 

0.66 

0.64 

0.32 

Ellrod-2 

4x10'^ 

0.61 

0.66 

0.63 

0.59 

0.29 

Endlich 

0.225 

0.63 

0.66 

0.53 

0.55 

0.19 

ITFA 

0.13 

0.57 

0.62 

0.67 

0.65 

0.29 

Richardson  # 

4 

0.58 

0.63 

0.66 

0.64 

0.29 

Essentially,  the  reports  were  broken  into  two  eategories:  one  eoneerned  with  any 
turbulenee  and  another  with  turbulenee  of  moderate  or  greater  (MOG)  intensity.  The 
results  are  listed  in  Table  2.  The  data  used  in  Table  2  was  adapted  from  Brown  et  al. 
(2000)  and  was  collected  between  21  December  1988  and  31  Mareh  1999  for  flights 
above  20,000  ft.  Table  2  shows  that  DTF3,  Ellrod-1  and  -2  and  ITEA  performed  the  best 
overall.  However,  it  is  apparent  the  difference  between  PODy  values  for  most  of  the 
indiees  used  in  this  study  was  relatively  small.  This  result  may  be  attributable  to  the 
strong  similarity  in  the  terms  involved  in  the  underlying  equations  for  the  various  indices. 
Additionally,  this  table  highlights  the  measurable  improvement  in  skill  every  index 
experiences  when  predicting  MOG  turbulence. 
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2.1.1.  Knapp-Ellrod  Index.  To  predict  turbulence  caused  by  Kelvin-Helmholtz 
instabilities,  Knapp  and  Ellrod  developed  two  indices  known  as  Ellrod-1  and  Ellrod-2 
(Ellrod  and  Knapp  1991).  These  indices  are  almost  identical;  however,  because  the 
convergence  term  tends  to  be  an  order  of  magnitude  smaller  than  the  other  terms,  Ellrod- 
1  scales  it  out.  The  indices  are  found  by  computing  the  product  of  the  sum  of 
deformation  and  convergence  at  a  pressure  level  with  vertical  shear  (Ellrod  and  Knapp 
1991).  The  indices  are  derived  from  a  version  of  Petterson’s  (1956)  equation  for 
frontogenic  intensity: 


4  =  0.5 


ydnj 


{DEF{cos2/3)  +  CVG) 


DEE  =  y!  DST^  +  DSH^ 


DST  = 


du  dv 
dx  dy 


dx  dy 


CVG  =  - 


du  dv 

- 1 - 

dx  dy  j 


(3) 

(4) 

(5) 

(6) 

(V) 


where  Aj.  is  the  time  rate  of  change  of  the  temperature  gradient  on  a  constant  pressure 
3  T 

surface,  —  is  the  magnitude  of  the  temperature  gradient  normal  to  the  isotherms,  DEE  is 
dn 

the  square  root  of  the  sum  of  the  squares  of  the  stretching  (DST)  and  shear  (DSH) 
deformation,  jd  is  the  angle  between  the  isotherms  and  the  axis  of  dilatation  (stretching 
axis  in  deformation  flow),  and  CVG  is  convergence.  The  temperature  gradient  term  is 
related  to  vertical  wind  shear  by  the  thermal  wind  equation: 
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where  f  is  the  eoriolis  parameter,  T  is  the  temperature,  g  is  the  gravitational  aeeeleration 


and 


ydZy 


is  the  vertieal  wind  shear.  After  retaining  the  signifieant  terms,  equation  (3) 


simplifies  to: 


77j  =  VWS(DEF) 

(9) 

TI^  =  VWS{DEF  +  CVG) 

(10) 

To  determine  the  thresholds  for  the  indiees,  Ellrod  and  Knapp  (1991)  eompared 
foreeasted  events  of  turbulence  to  PIREPS.  Table  3  lists  the  resulting  thresholds.  Air 
Eorce  Global  Weather  Central  (AEGWC),  now  known  as  the  Air  Eorce  Weather  Agency 
(AEWA),  calculated  TI^ ,  or  Ellrod-2,  from  the  AEGWC  global  spectral  model  (GSM) 
and  High-Resolution  Analysis  System  (HIRAS).  The  National  Meteorological  Center 
(NMC)  calculated  TYj ,  or  Ellrod- 1,  from  the  National  Meteorological  Center  Nested  Grid 

Model  and  Aviation  Model  to  determine  their  values. 

Table  3,  Turbulence  intensity  threshold  for  Ellrod  index 


Turbulence  Intensity 

AFGWC 

NMC  NGM 

NMC  AVN 

LGT-MDT 

4 

MDT 

8 

4 

2 

MDT-SVR 

12 

8 

4 

The  contingency  table  (see  Eigure  1)  verification  statistics  at  both  AEGWC  and 
NMC  showed  some  statistical  skill  for  several  criteria.  Eor  this  study,  the  probability  of 
detection  (POD)  is  the  PC  as  outlined  in  Table  1.  The  POD  was  above  70%  for  both 


AFGWC  and  NMC.  The  FAR  (as  described  in  Table  1)  was  22%  for  the  aviation  model 
(AVN),  36%  for  the  nested  grid  model  (NGM)  and  30-48%  for  AFGWC. 

It  is  important  to  remember  the  mesoscale,  short-lived  nature  of  CAT  and  the 
relatively  poor  pilot  report  sampling  of  CAT  identified  regions  when  considering  the  high 
probability  of  false  alarms.  These  indices  predict  regions  where  conditions  are  favorable 
for  CAT  to  develop  and  not  where  it  is  actually  occurring.  The  problem  with  PIREP 
sampling  is  discussed  further  in  Section  2.3. 

Additionally,  in  a  comparative  study  from  November  1989  to  April  1990,  the 
Ellrod  Index  compared  well  with  the  Empirical  Dutton  Index,  both  having  similar  POD 
and  FAR  (Ellrod  and  Knapp  1991).  The  Dutton  index  was  determined  using  linear 
regression  software  from  the  University  of  California  Los  Angeles  (Dutton  1980). 

Dutton  compared  PIREP  data  to  several  potential  indicators  of  CAT  including  the  rate  of 
change  of  Richardson  number  (following  the  flow),  horizontal  wind  shear,  vertical  wind 
shear,  wind  speed,  vertical  velocity,  horizontal  gradient  of  vertical  velocity,  vorticity,  and 
deformation.  Using  the  regression  software,  it  determined  an  equation  using  the  sum  of 
the  square  of  vertical  wind  shear  and  horizontal  wind  shear  as  the  best  predictor  of 
turbulence. 

2.2  Mountain  Waves 

Mountain  waves  have  long  been  recognized  as  a  cause  for  CAT.  They  typically 
form  when  the  normal  component  of  the  wind  to  the  mountain  is  10  m/s  or  more.  They 
are  known  to  propagate  200km  downstream  of  the  mountain  that  causes  them  and  their 
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intensity  is  usually  maximized  1300m  above  and  below  the  tropopause  (Hopkins  1977). 
Figure  2  shows  a  elassic  model  of  a  mountain  wave  and  areas  of  turbulenee  produeed  by 
it  (Chandler  1987). 


Mountain  waves  are  a  type  of  internal  gravity  wave.  An  internal  gravity  wave  is 
caused  by  the  combined  force  of  gravity  and  buoyancy  on  a  parcel  of  air  displaced 
vertically  (NRL  2003).  When  displaced,  the  forces  of  gravity  and  buoyancy  are  typically 
no  longer  in  balance.  This  imbalance  causes  the  parcel  to  accelerate  vertically  until  it 
reaches  a  level  where  the  forces  are  again  in  balance.  In  an  unstable  environment,  an 
upward  displaced  air  parcel  is  warmer  than  its  surroundings.  This  parcel  will  continue  to 
accelerate  upward  away  from  the  original  layer.  In  a  stable  environment,  an  upward 
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displaced  air  parcel  is  eooler  than  its  surroundings.  In  this  ease  the  foree  of  gravity  is 
greater  than  buoyaney  and  the  parcel  aecelerates  downward  toward  (and  beyond)  its 
original  level.  This  aeeeleration  past  its  original  level  leaves  it  warmer  than  its 
surroundings  and  now  the  buoyancy  force  overcomes  gravity.  This  causes  the  parcel  to 
aeeelerate  upward  towards  its  original  level.  The  situation  repeats  itself  causing  the 
pareel  to  oscillate  vertieally  and  laterally  about  some  point  in  the  horizontal  plane  (NRL 
2003).  In  the  case  of  mountain  waves,  this  motion  is  continually  forced  by  a  mountain 


range. 

The  frequeney  with  which  the  air  pareel  oseillates  vertieally  is  known  as  the 
Brunt-Vaisala  frequeney  and  is  expressed  mathematically  as: 


where  N  is  the  Brunt-Vaisala  frequeney  and  6  is  the  potential  temperature.  The 
magnitude  of  the  Brunt-Vaisala  frequency  is  based  on  the  statie  stability  of  the 
environment.  N  also  makes  up  the  numerator  in  computing  the  Richardson  number  for 
Kelvin-Helmholtz  instabilities. 


Based  on  the  discussion  in  Holton  (1992),  the  dispersion  relation  relates  the 
frequency  of  oscillation  to  the  wave  number.  This  relationship  determines  whether  the 
phase  speed  is  a  funetion  of  wavenumber  and  whether  the  wave  pattern  ehanges  shape 
(disperses).  In  the  ease  of  topographieally  foreed  gravity  waves  whieh  are  stationary  with 
respeet  to  the  ground,  the  dispersion  relation  is: 

m  =—^-k  (12) 

u 
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where  m  is  the  vertieal  wavenumber,  u  is  the  mean  wind  flow  perpendieular  to  the 
mountain  range,  and  k  is  the  horizontal  wavenumber. 

In  equation  12,  it  is  possible  for  m  to  be  less  than  zero.  In  this  ease,  the  waves 
are  vertieally  trapped.  Typieally,  flow  over  narrow  ridges  produees  vertieally  trapped 
waves.  Narrow  ridges  tend  to  cause  smaller  wavelengths  and  thus  larger  values  of 
horizontal  wavenumber.  Additionally,  small  stability  or  a  large  perpendicular  component 
of  the  wind  might  cause  vertically  trapped  waves.  In  either  case,  vertically  trapped  waves 
typically  decay  exponentially  with  height. 

When  m  is  greater  than  or  equal  to  zero,  the  waves  are  vertically  propagating. 
Vertically  propagating  waves  occur  with  some  combination  of  small  mean  perpendicular 
flow,  large  stability,  and  large  horizontal  wavelength.  Vertically  propagating  waves  must 
transport  energy  upward  because  the  ground  prevents  downward  propagation.  Vertically 
trapped  waves  also  propagate  energy  upward  but  the  energy  decays  with  height.  Also,  it 
is  possible  for  both  wave  types  to  exist  simultaneously  (Neiman  and  Shaw  2003).  In  any 
case,  wave  energy  propagates  upward  away  from  the  ground,  although  vertically  trapped 
wave  amplitude  decreases. 

In  the  vertically  propagating  case,  the  amplitude  of  the  wave  grows  with  height  as 
the  air  becomes  less  dense  (NRL  2003).  Occasionally,  the  wave  breaks  and  turbulence 
results.  The  severity  of  the  turbulence  depends  on  the  amplitude  of  the  wave  when  it 
breaks.  The  amplitude  growth  with  height  explains  why  mountain  wave  turbulence  tends 
to  become  more  severe  at  higher  altitudes. 

Until  recently,  forecasters  have  relied  on  rules  of  thumb  and  rawinsonde 
observation  (RAOB)  data  to  predict  the  occurrence  of  mountain  waves.  When  using 
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RAOB  data,  they  assessed  the  likelihood  of  turbulenee  based  on  software  sueh  as  the 
HiCAT  (High  Altitude  CAT)  analysis  module  whieh  looks  for  “S”  shaped  temperature 
patterns  in  three  adjaeent,  distinet  layers  and  eomputes  the  parameter  based  on  the  lapse 
rate,  mixing  layer,  depth  of  the  “S”  layer  and  the  “S”  layer  vertieal  temperature  differenee 
(Sinelair  and  Kuhn  1991).  Reeently,  numerieal  model  algorithms  sueh  as  the  Mountain 
Wave  Foreeast  Model  (MWFM)  have  been  developed  and  show  potential  for  predieting 
mountain  wave  turbulenee  in  the  troposphere  (Eekerman  et  al.  2000). 

2.2.1  Mountain  Wave  Forecast  Model  VI. 1.  Dr.  Baemeister  (1993)  of  the  Naval 
Research  Laboratory  introduced  MWFM  VI .  1 .  MWFM  VI .  1  calculates  the  distribution 
of  mountain  waves  forced  by  air  flow  over  a  set  of  ridges  contained  in  a  ridge  database. 
The  ridge  database  was  generated  using  a  global  topography  dataset,  compiled  by  the 
National  Center  for  Atmospheric  Research  (NCAR).  This  dataset  was  used  to  determine 
ridge  latitude,  longitude,  height,  width,  and  orientation.  MWFM  Vl.l  combines  this 
information  and  a  two-dimensional  hydrostatic  gravity  wave  model  to  determine  the 
severity  and  extent  of  the  waves  generated  by  each  ridge.  For  a  detailed  description  of 
the  algorithm  used  in  determining  the  ridge  database,  refer  to  Baemeister  et  al.  (1993). 

According  to  Baemeister  et  al.  (1994),  MWFM  Vl.l  derives  a  turbulence  estimate 
using  the  following  method.  First,  the  following  assumptions  are  made:  wave  crests 
generated  by  a  ridge  remain  parallel  to  the  ridge,  the  waves  generated  by  the  flow  over 
the  ridge  are  hydrostatic,  and  the  average  momentum  flux  of  the  waves  remains  constant 
with  height  until  the  wave  breaks.  The  first  calculation  uses  the  perpendicular  component 
of  the  wind  to  each  ridge  to  determine  the  stratification  frequency  profile  over  that  ridge. 
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Next,  the  first  guess  wave  displaeement  at  eaeh  level  is  ealeulated  using  the  following 
equation: 


dl{z  +  ISz)  =  d^{z) 


p{z)N k{z)U  ^.,{z) 
p{z  +  Az)Ni^  (z  +  Az)U^.^  (z  +  Az) 


(13) 


whieh  is  basieally  a  ratio  of  the  produet  of  density  p ,  stratifieation  frequeney  N^. ,  and 
perpendieular  eomponent  of  the  wind  U at  two  adjaeent  levels.  If  this  value  exeeeds 

the  saturation  limit,  it  is  redueed  to  equal  the  saturation  limit.  To  eompute  the  wave 
displaeement  at  the  surfaee,  the  minimum  of  the  saturation  limit  at  the  surfaee  and  a 
multiple  of  the  ridge  height  parameter  is  determined.  The  ridge  height  parameter  is  based 
on  the  deviation  of  the  ridge  height  from  the  surrounding  topography.  That  value  is 
multiplied  by  a  faetor  of  four  for  the  eomparison. 

Whenever  the  wave  amplitude  approaehes  or  exeeeds  the  saturation  limit,  some  of 
the  energy  is  assumed  to  eontribute  to  eonveetive  and  shear  instability.  This  results  in 
turbulenee.  Turbulenee  intensity  is  the  differenee  between  momentum  flux  eomputed  for 
two  adjaeent  levels.  Momentum  flux  is  proportional  to  the  terms  in  equation  13.  Due  to 
the  relatively  simple,  inexpensive  ealeulations  required,  MWFM  VI.  1  is  able  to  produee 
timely,  world-wide  foreeasts  for  altitudes  up  to  the  mesosphere  (Eekerman  et  al.  2000). 
For  a  more  detailed  look  at  the  equations  involved  in  ealeulating  MWFM  Vl.l,  refer  to 
Appendix  A. 


2.2.2.  Mountain  Wave  Forecast  Model  V2.1.  MWFM  V2.1  addressed  some  of  the 
simplifieations  made  by  V 1 . 1 ,  at  the  expense  of  eomputation  time,  in  an  effort  to  more 
aeeurately  model  mountain  wave  dynamies.  The  primary  ehanges  ineluded  the 
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conversion  from  two  to  three  dimensions  and  from  hydrostatie  to  non-hydrostatie 
caleulations.  Although,  the  model  was  primarily  designed  for  stratospherie  turbulenee,  it 
has  shown  promise  for  use  in  the  middle  and  upper  portions  of  the  troposphere 
(Eekerman  et  al.  2000). 

Aeeording  to  Marks  and  Eekermann  (1995),  MWEM  V2.1  is  ealeulated  using  a 
ray  traeing  model.  The  ray  traeing  model  tracks  the  trajeetory  and  amplitude  deviations 
of  individual  waves  propagating  through  the  atmosphere.  In  MWEM  V2. 1 ,  a  user- 
defined  number  of  rays  are  launehed  from  eaeh  ridge  with  a  pattern  that  depends  on  the 
ratio  of  ridge  width  to  length  (Eekermann  et  al.  2001).  In  the  ease  of  a  wide  ridge,  the 
rays  are  launehed  at  equidistant  angles  from  0-180°  on  the  downwind  side  of  the  ridge. 
The  eloser  a  ray  is  to  the  90°  angle,  the  higher  the  initial  magnitude  it  is  given.  In  the 
ease  of  an  isolated  peak,  the  rays  take  on  a  pattern  similar  to  one  ereated  by  a  ship’s  wake 
and  the  rays  are  assigned  similar  initial  amplitudes.  Typieally,  18  rays  are  launehed  from 
eaeh  ridge  with  the  MWEM  V2. 1 . 

Eaeh  ray  is  assigned  initial  values  of  longitude,  latitude,  and  altitude. 

Additionally,  the  horizontal  wavenumber,  ground-based  frequency  and  other  terms  are 
ehosen  sueh  that  the  vertieal  group  veloeity  is  positive  (Marks  and  Eekermann  1995). 
This  ensures  the  wave  energy  is  propagated  upward.  Next,  the  wave  amplitude  is 
determined  by  ealeulating  the  vertieal  flux  of  wave  aetion.  This  flux  is  determined  by 
multiplying  the  vertieal  group  veloeity  and  the  wave  aetion  density  together.  The  wave 
aetion  density  is  simply  the  wave  energy  density  divided  by  the  intrinsie  (Eagrangian) 
frequeney  of  the  wave.  The  vertical  flux  of  wave  aetion  is  then  integrated  to  find  the 
ehange  in  the  vertieal  flux  of  wave  aetion,  whieh  is  used  to  find  the  wave  amplitude 
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(Allen  2003).  Once  the  wave  amplitude  is  known,  the  calculations  discussed  in  section 
2.2.1  are  applied  to  determine  the  turbulence  intensity.  For  a  detailed  description  of 
MWFM  V2.1,  refer  to  Marks  and  Eckermann  (1995). 

According  to  Eckermann  as  cited  in  Allen  (2003),  the  MWEM  V2.1  also  accounts 
for  the  effects  of  friction  and  blocking  using  an  arbitrary  scaling  factor  and  the  local 
Eroude  number  respectively.  An  arbitrary  scaling  factor  accounts  for  friction  by  reducing 
the  wind  speed  in  the  calculations.  The  local  Eroude  number  is  the  ratio  of  kinetic  energy 
to  potential  energy.  When  the  Eroude  number  is  on  the  order  of  one,  the  effective  height 
of  the  terrain  and  the  wave  amplitude  are  reduced  to  account  for  blocking  effects. 

2.3  Pilot  Reports 

PIREPS  are  commonly  considered  the  best  method  available  for  verifying 
turbulence  forecasts.  These  reports  are  typically  reported  by  the  aircrew  over  the  radio, 
but  some  aircraft  have  equipment  which  automatically  records  weather  information. 
Aircrews  typically  base  the  intensity  of  the  turbulence  they  report  on  their  subjective 
interpretation  while  automated  reports  typically  determine  turbulence  intensity  by 
measuring  the  center  of  gravity  acceleration  peaks.  Different  sources  have  related 
different  values  of  these  peaks  to  different  turbulence  intensity.  Current  World 
Meteorological  Organization  (WMO)  standards  use  the  following  thresholds:  light  (0.15 
g  <  g*  <0.5  g),  moderate  (0.5  g  <  g*  <  l.Og),  and  severe  (g*  >  1.0  g)  to  quantify  CAT 
(WMO-No.  958  2003).  These  values  are  based  on  g*  being  the  acceleration  deviations 
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from  the  standard  1 .0  g  acceleration.  The  system  reports  turbulence  measurements  from 
these  sensors  as  turbulence  type  while  aircrews  report  turbulence  intensity. 

PIREPS  of  turbulence  (or  lack  of  turbulence)  are  typically  entered  into  a  2X2 
contingency  table  with  forecasts  of  turbulence  (or  lack  of  turbulence)  as  in  Figure  1 .  The 
values  in  these  tables  are  then  used  with  equations  such  as  those  listed  in  Table  1  to 
summarize  the  skill  of  a  forecast  method  (Brown  and  Young  2000). 

Despite  being  the  best  verification  tool  available,  using  PIREPS  does  pose  several 
important  problems.  First,  these  reports  frequently  rely  on  a  pilot’s  subjective 
interpretation  of  the  severity  of  the  turbulence  encountered  instead  of  relying  on  an 
objective,  quantitative  value.  This  may  lead  to  inconsistencies  in  the  reporting  of  CAT 
severity.  Second,  aircraft  do  not  systematically  and  completely  sample  a  grid  to 
definitively  determine  the  presence  of  CAT  (Brown  and  Young  2000).  Essentially,  an 
aircraft  usually  flies  a  fairly  direct  route  between  two  points  instead  of  flying  a  grid 
pattern  over  a  region.  A  simulation  study  conducted  by  Brown  and  Young  (2000) 
demonstrates  the  inaccuracy  in  computing  quantities  such  as  FAR,  bias,  and  other 
statistics  when  using  PIREPS  for  CAT  verification.  In  that  study,  various  percentages 
of  “yes”,  “no”,  or  “yes”  and  “no”  PIREPS  were  removed  from  a  sample  of  10,000 
PIREPS  and  the  various  indices  were  recomputed.  As  larger  percentages  of  PIREPS 
were  removed,  the  PODy  and  PODn  values  remained  relatively  constant.  However,  the 
FAR,  bias,  and  CSI  changed  significantly  when  a  portion  of  either  the  “yes”  or  “no” 
PIREPS  were  removed  (not  both).  Removal  of  “yes”  PIREPS  resulted  in  increased  FAR, 
bias,  and  CSI  while  removal  of  “no”  PIREPS  resulted  in  decreased  FAR,  bias,  and  CSI. 
This  suggests  a  strong  dependence  on  the  relative  frequency  of  “yes”  or  “no”  PIREPS  to 
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the  FAR  (resulting  from  PIREP  sample  having  different  probability  distribution  function 
than  the  actual  turbulence)  and  that  high  FAR  is  less  significant  than  it  might  initially 
seem. 

There  are  several  potential  problems  with  accelerometer  data  although  they  are 
more  objective  than  PIREPS.  First,  depending  on  the  size,  altitude,  and  airspeed  of  the 
aircraft  or  the  nature  of  the  turbulence,  the  accelerometer  may  read  different  values  on  the 
accelerometer  for  the  same  turbulence  encounter.  Furthermore,  the  pilot  may  take  actions 
that  modify  the  attitude  of  the  aircraft  and  cause  the  accelerometer  value  to  change. 

These  two  effects  limit  the  application  of  AVAR  data  as  a  tool  in  determining  turbulence 
intensity. 
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III.  Methodology 


3.1  Overview 

The  primary  goal  of  this  research  was  to  provide  the  weather  forecaster  with  a 
quick,  accurate  method  of  developing  turbulence  charts  that  incorporate  both  primary 
causes  of  CAT.  To  best  achieve  that  goal,  the  needs  for  accuracy  and  speed  were 
addressed  separately.  This  provided  one  distinct  advantage.  The  speed  with  which  the 
turbulence  charts  get  developed  is  independent  of  the  quality  of  the  underlying  indices. 
Clearly,  the  forecaster  aids  must  be  as  skillful  as  possible  to  produce  the  most  accurate 
charts.  Over  time,  the  quality  and  resolution  of  the  model  might  improve  or  a  superior 
new  index  might  be  developed.  Regardless,  the  actual  method  of  developing  the  charts 
does  not  have  to  change,  making  implementation  of  improvements  seamless  to  the 
forecaster  building  the  charts. 

To  capture  turbulence  caused  by  Kelvin-Helmholtz  instabilities,  the  Ellrod-2 
index  was  chosen.  This  index  was  chosen  for  two  reasons.  First,  past  research  has 
identified  the  Ellrod-2  index  as  being  relatively  skillful  in  predicting  turbulence  (Brown 
et  al.  2000).  Second,  the  Ellrod-2  index  is  built  into  the  software  currently  used  by  USAF 
weather  squadrons.  This  allows  for  easy  integration  of  this  method  into  OWS  operations. 

The  Mountain  Wave  Forecast  Model  was  chosen  to  identify  turbulence  caused  by 
mountain  waves.  As  discussed  previously,  there  are  two  versions  of  the  model  and  for 
this  study  both  versions  were  tested.  The  MWFM  was  chosen  for  two  reasons  also. 
Primarily,  this  model  has  shown  some  potential  in  predicting  turbulence  in  the 
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troposphere  caused  by  mountain  waves  (Eckermann  et  al.  2000).  Additionally,  the 
MWFM  is  already  being  implemented  at  AFWA  to  predict  stratospheric  turbulence.  The 
modifications  required  to  implement  the  model  at  tropospheric  levels  should  be  relatively 
easy. 

The  Ellrod-2  index  and  MWFM  indices  were  calculated  in  a  window  from  20N- 
60N  and  140W-30W.  All  but  29  of  the  analyses  from  May  20,  2003  to  October  15,  2003 
were  used,  along  with  PIREPS,  to  determine  the  thresholds.  Additionally,  the  12-hour 
and  24-hour  forecast  data  were  used  to  make  forecasts  from  October  21,  2003  to  October 
27,  2003  and  from  November  5,  2003  to  November  12,  2003.  These  forecasts  were 
verified  against  PIREPS.  Finally,  12-hour  and  24-hour  operational  forecasts  produced  by 
USAF  weather  squadrons  from  October  21,  2003  to  October  27,  2003  and  from 
November  5,  2003  to  November  12,  2003  were  also  verified  with  PIREPS. 

3.2  Data 

3.2.1  Model  Input  Data.  The  GFS  model  data  were  downloaded  from  the  NCEP  FTP 
site.  The  GFS  model  is  run  four  times  daily  although  this  research  only  used  the  OOZ  and 
12Z  model  runs.  The  model  also  forecasts  out  to  120  hours.  This  research  only  looked  at 
the  analysis,  12-hour  forecast,  and  24-hour  forecast.  Horizontally,  the  GFS  model  uses 
spectral  basis  functions  with  resolution  T254.  This  amounts  to  approximately  0.5°  x  0.5° 
grid  point  spacing  in  the  horizontal.  Vertically,  the  GFS  model  is  composed  of  64  sigma 
layers  with  emphasis  below  800  mb  and  above  100  mb.  This  research  used  the  layers 
between  800  mb  and  100  mb.  The  data  files  used  in  this  study  consisted  of  l°xl°  global. 
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gridded  model  data  in  the  gridded  binary  (GRIB)  format.  More  detail  on  the  GFS  model 
is  available  at  the  NCEP  web  site  (NCEP  2003).  The  ETP  was  performed  by  either  the 
MWEM  software  or  by  other  methods  that  stored  the  data  locally.  The  MWEM  was 
called  from  a  script  which  first  determined  if  the  data  were  available  locally  before 
attempting  to  retrieve  the  data  from  NCEP  using  ETP. 

The  retrieved  data  were  stored  on  a  36  gigabyte  hard  drive.  To  reduce  the  storage 
space  required  by  these  large  data  fdes,  the  data  was  compressed  by  GNU  zip  software 
available  for  download  on-line.  Compression  reduced  the  file  sizes  to  about  80%  of  their 
original  size.  This  resulted  in  a  saving  of  about  5  megabytes  (MB)  on  an  approximately 
26  MB  file  retrieved  from  AEWA.  The  data  stored  in  a  local  archive  was  retrieved  from 
the  same  ETP  site  but  the  files  were  trimmed  to  include  only  0-90N  latitude  and  180W-0 
longitude  instead  of  the  entire  earth.  Compression  on  these  files  resulted  in  a  savings  of 
about  1  MB  on  a  5.5  MB  file.  Using  the  5.5  MB  data  files  required  modification  of  the 
MWEM  and  Ellrod-2  index  code  to  account  for  the  different  starting  and  ending  latitude 
and  longitude  associated  with  those  files.  Einally,  a  third  set  of  data  files  only  included 
the  necessary  model  fields  for  the  necessary  levels  on  the  global  grid.  These  files  were 
reduced  from  about  10  MB  to  about  8  MB  using  compression  software. 

3.2.2  MWFM Input  Data.  The  MWEM  contains  a  variable  that  determines  the  source  of 
the  model  data.  Based  on  data  availability,  the  “ftpoption”  variable  was  used  with  either 
a  “0”  option  or  a  “2”  option.  The  “0”  option  causes  the  model  to  retrieve  its  data  locally 
and  not  to  erase  the  data  after  the  model  run  is  complete.  The  “2”  option  causes  the 
model  to  retrieve  the  data  from  the  NCEP  ETP  site  and  not  to  erase  the  data  after  the 
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model  run  is  complete.  The  MWFM  used  four  fields  from  the  GFS  model.  These 
include  the  absolute  temperature,  geopotential  height,  U  wind  component,  and  V  wind 
component.  These  fields  were  used  by  the  MWFM  equations  to  calculate  turbulence  for 
every  50  millibar  (mb)  layer  between  100  mb  and  750  mb  with  the  highest  layer  being  the 
100-150  mb  layer  and  the  lowest  being  the  700-750  mb  layer.  These  layers  encompass 
the  entire  section  of  the  atmosphere  necessary  to  produce  turbulence  charts  from  FLIOO 
to  FL500. 


3.2.3  Ellrod-2  Input  Data.  The  Ellrod-2  index  was  calculated  using  the  same  GFS  model 
files  used  by  the  MWFM.  The  GRIB  files  were  converted  to  binary  files  by  a  UNIX  shell 
script  to  allow  a  formula  translation  (FORTRAN)  program  to  read  them.  The  shell  script 
extracted  the  U  wind  component,  V  wind  component,  and  geopotential  height  every  50 
mb  from  50-800  mb  and  saved  each  as  an  individual  binary  file.  These  binary  files  were 
used  to  calculate  the  Ellrod-2  index  from  the  750  mb  level  to  the  150  mb  level  and  to 
determine  the  average  height  over  the  same  layers  as  the  MWEM. 

3.2.4  PIREP  Input  Data.  A  large  database  of  military  and  civilian  PIREPS  covering  the 
available  data  was  received  from  the  Air  Eorce  Combat  Climatology  Center  (AECCC). 
This  database  included  the  latitude,  longitude,  year,  month,  day,  hour,  altitude,  and 
turbulence  intensity  fields.  Unfortunately,  the  category  (or  type)  of  aircraft  was  not 
available  in  this  database.  As  a  result,  all  pilot  report  turbulence  intensities  were  treated 
as  if  they  originated  from  a  category  2  aircraft,  since  weather  squadron  turbulence 
forecasts  are  geared  towards  that  category  of  aircraft.  Eor  each  model  run  time,  any 
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PIREPS  within  a  6  hour  window  centered  on  the  valid  time  was  saved  in  a  text  file.  This 


text  file  was  passed  to  the  FORTRAN  program  for  further  analysis. 

This  database  represented  turbulenee  intensity  with  a  number.  Aecording  to 
AFMAN  15-124  (2001)  and  WMO-No.  306  (1995)  code  table  0300,  the  number  0 
eorresponds  to  none,  1  to  light,  2-5  to  moderate,  6-9  to  severe,  and  “X”  to  extreme. 
Occasional  turbulence  suggests  turbulence  oceurring  less  than  1/3  of  the  time.  Table  4 
details  the  turbulenee  intensities  as  outlined  in  AFMAN  15-124  and  WMO  standards. 
These  intensities  were  later  sealed  so  that  values  of  0  were  assigned  0,  1  were  assigned  1, 
2-5  were  assigned  2,  and  6-9  were  assigned  3.  All  told,  a  sample  of  18,165  PIRFPS  were 
used  in  this  researeh  study,  of  whieh  15,745  (May  20,  2003  -  Oetober  15,  2003)  were 
used  to  set  the  thresholds  and  2,420  (October  21,  2003  -  Oetober  27,  2003,  November  5, 
2003  -  November  12,  2003)  were  used  to  verify  the  forecasts. 


Table  4.  List  of  turbulence  intensity  values 


Value 

Turbulence  Intensity 

Value 

Turbulence  Intensity 

0 

None 

6 

Occasional  severe 

1 

Fight 

7 

Frequent  severe 

2 

Occasional  moderate 

8 

Oecasional  severe  in  eloud 

3 

Frequent  moderate 

9 

Frequent  severe  in  eloud 

4 

Occasional  moderate  in  eloud 

X 

Extreme 

5 

Frequent  moderate  in  cloud 

5.5  Index  Configuration 


3.3.1  MWFM  Configuration.  For  this  experiment,  both  versions  of  the  MWFM  were 
tested.  For  both  versions,  the  MWFM  program  suite  is  set  up  in  modules  making  it  easy 
to  modify  experiments  for  different  models,  locations  on  the  earth,  and  regions  of  the 
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atmosphere.  In  this  study  the  model  was  run  over  the  entire  eontinental  United  States 
(CONUS)  from  a  window  ineluding  20N-60N  latitude  and  140W-30W  longitude.  This 
window  provided  diverse  topography  and  a  variety  of  weather  eonditions  to  test  the 
indie  es. 

This  study  initiated  the  MWFM  from  a  shell  seript  whieh  passed  the  neeessary 
date  information  to  start  the  program,  moved  the  loeal  arehive  files  to  the  appropriate 
MWFM  data  direetories,  ehose  whether  to  run  “ftpoption”  0  or  2,  moved  the  output  to  the 
appropriate  output  direetory,  and  ensured  all  the  neeessary  files  were  stored  in  the  proper 
loeation.  When  initiated,  the  MWFM  starts  the  Interaetive  Data  Language  (IDL) 
interface.  It  then  immediately  calls  the  at  glider.pro  IDL  script.  This  script  controls  the 
forecasting  sequence  by  calling  any  necessary  IDL  scripts  or  FORTRAN  programs 
depending  on  several  input  parameters.  Most  notably,  the  at  glider.pro  script  calls  the 
what2do  file  which  outlines  the  experiment  and  drv  glider.pro  which  calls  the  actual 
forecast  routines. 

The  what2do  file  determines  where  and  how  the  MWFM  will  run  an  experiment. 

It  also  outlines  how  any  output  will  look  and  what  parameters  the  program  should 
provide  to  the  user.  It  is  one  of  several  files  that  must  be  manipulated  to  set  up  a  new 
experiment.  For  this  experiment,  the  what2do  fide,  what2do_belsontest.pro,  was  set  up 
to  run  the  model  over  the  continental  United  States  (CONUS)  for  all  pressure  levels  that 
comprise  the  10,000  ft  to  50,000  ft  MSL  altitude  range.  The  levels  and  field  keyword 
were  set  up  to  reflect  the  desire  to  view  the  turbulence  intensity  variable  for  all  the  50 
millibar  layers  between  the  800-850  mb  layer  and  the  100-150  mb  layer.  The  latitude  and 
longitude  limits  variable  was  set  up  to  cover  140W  to  SOW  longitude  and  20N  to  60N 
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latitude.  Furthermore,  the  what2do  file  was  set  up  to  run  both  VI .  1  and  V2. 1  with  GFS 
model  data  as  input.  For  a  copy  of  the  what2do  file  used  in  this  experiment,  refer  to 


Appendix  B. 


Turbulence  Intensities  (EP  Flux  Deposition)  ^(z)-<)t(z+Az) 
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Figure  3.  Example  of  output  from  MWFM  V2. 1 


The  output  from  MWFM  VI .  1  and  V2. 1  was  set  to  save  graphical  and  raw  data.  The 
graphical  output  was  saved  in  postscript  format.  It  was  only  used  to  test  the  accuracy  of 
the  FORTRAN  code  which  converts  the  MWFM  output  to  gridded  output  for  display 
with  the  Grid  Analysis  and  Display  System  (GrADS)  software  and  was  not  configured 
for  ideal  viewing.  An  example  of  the  graphical  output  used  for  testing  is  provided  in 
Figure  3.  In  Figure  3,  the  dark  region  along  the  Califomia/Nevada  border  indicates 
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turbulence.  Additionally,  two  scripts  saved  the  raw  data  with  latitude,  longitude,  and 
turbulence  intensity  for  each  ridge  in  the  window. 

3.3.2  Ellrod-2  Index  Configuration.  The  Ellrod-2  index  was  implemented  with  the  same 
equations  used  in  calculating  the  index  at  the  Air  Force  Weather  Agency.  The  equations 
used  2"‘*  order  centered  finite  differencing  to  compute  the  vertical  wind  shear  (VWS), 
convergence  (CVG),  and  deformation  (DEF)  terms.  Elsing  centered  differencing  when 
computing  VWS  introduces  some  error  since  the  difference  in  height  between  pressure 
levels  is  not  constant  and  the  center  pressure  level  is  not  in  the  geometric  center.  This 
simplification  resulted  in  an  error  of  5%  or  less  in  the  computation.  Also,  geopotential 
height  was  substituted  for  geometric  height  in  the  computations.  According  to  Holton 
(1992),  the  geometric  height  and  geopotential  height  are  almost  identical  in  the  regions  of 
interest  for  this  study.  Computing  the  index  for  a  level  requires  the  U  and  V  wind  data 
for  three  adjacent  levels,  centered  on  the  desired  level.  It  also  required  the  geopotential 
height  for  the  levels  directly  above  and  below  the  desired  level.  Distances  were 
calculated  using  the  haversine  method,  as  described  by  Sinnott  (1984).  The  equations 
used  in  the  haversine  method  are  provided  in  Appendix  A  and  the  FORTRAN  subroutine 
is  provided  in  Appendix  E,  subroutine  Haversine.  For  the  FORTRAN  code  used  to 
compute  the  Ellrod-2  index,  refer  to  Appendix  E  subroutine  Ellrod2. 
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3.4  Implementation 


The  indices  were  combined  using  a  series  of  shell  scripts  and  a  FORTRAN 
program.  The  controlling  script  retrieved  the  necessary  analysis,  12-hour  forecast,  and 
24-hour  forecast  data  files  (GFS  U,  V,  and  GPH  fields,  MWFM  Vl.l  and  V2.1  raw  data 
files,  and  PIREPS  based  on  the  date)  required  for  the  FORTRAN  program.  For  each 
level,  the  FORTRAN  program  computed  the  Ellrod-2  index  with  the  GFS  data  fields,  the 
MWFM  with  the  raw  data,  the  average  layer  height,  and  the  index  value  for  each  pilot 
report  with  the  truncated  PIREPS  text  file.  The  indices  were  saved  to  a  data  file  for 
display  using  GrADS  software  while  the  PIREPS  data  was  saved  to  a  text  file  for 
statistical  study. 

The  FORTRAN  program  also  created  another  chart  that  provided  a  composite 
view  of  the  indices.  This  chart  was  designed  to  provide  the  forecaster  with  a  snapshot  of 
the  atmosphere  from  above  that  would  allow  them  to  quickly  determine  the  horizontal 
location  of  the  turbulent  regions.  This  information  was  stored  in  the  same  data  file  as  the 
other  index  files. 

3.4.1  Index  Implementation.  The  Ellrod-2  index  was  computed  at  every  available  level 
between  100  mb  and  750  mb.  This  amounted  to  14  levels  separated  by  50mb  each. 

Since  this  computed  the  index  for  a  specific  level  and  the  MWFM  computed  the  index  for 
a  50  mb  layer,  it  was  impossible  to  make  the  indices  represent  the  exact  same  region  of 
the  atmosphere.  This  problem  was  addressed  by  assigning  the  bottom  of  the  MWFM 
layer  and  the  corresponding  Ellrod-2  index  to  the  same  chart.  For  example,  the  500-550 
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mb  MWFM  layer  and  the  550  mb  Ellrod-2  level  were  displayed  on  the  same  chart.  The 
choice  was  made  to  display  the  Ellrod-2  index  from  150  mb  to  750  mb  instead  of  from 
100  mb  to  700  mb  since  the  winds  at  the  top  of  the  column  are  above  the  troposphere  and 
typically  do  not  vary  as  much  as  the  winds  in  the  troposphere. 

Additionally,  the  average  geopotential  height  for  each  layer  was  computed.  This 
computation  was  performed  for  each  grid  point  in  the  model  by  averaging  the  height  for 
two  adjacent  levels.  The  results  were  used  as  the  heights  to  assign  values  of  each  index 
to  the  PIREPS  for  the  MWEM  indices.  The  Ellrod-2  index  levels  were  determined  using 
the  height  of  the  actual  level  for  which  the  Ellrod-2  index  was  computed.  The  PIREPS 
assignment  scheme  is  discussed  further  in  Section  3.4.2. 


Eigure  4.  Example  of  accuracy  of  MWEM  display  using  GrADS 

The  MWEM  indices  were  computed  for  each  grid  point  in  the  region  of  interest. 
The  raw  data  only  contained  a  subset  of  latitudes  and  longitudes  for  the  region  that 
matched  ridges  as  contained  in  the  MWFM  database.  This  raw  data  included  turbulence 
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intensities  greater  than  0.0.  A  representative  graphic  on  the  grid  was  made  using  the 
largest  value  of  the  index  in  a  1°  X  1°  box  centered  on  each  grid  point.  This  ensured  that 
the  maximum  value  was  assigned  to  every  point  in  the  array,  guaranteeing  the  worst-case 
is  always  chosen  (a  desirable  trait  where  flight  safety  is  concerned).  If  there  were  no 
values  of  raw  data  within  that  box,  the  grid  point  was  assigned  a  value  of  zero.  This 
method  provided  accurate  results  when  compared  to  the  postscript  fdes  produced  by  the 
MWFM.  Figure  4  provides  an  example  of  the  comparative  output  from  the  MWFM  and 
the  FORTRAN  code  with  MWFM  output  on  the  left.  The  peak  magnitudes  for  the  two 
images  in  Figure  4  are  located  in  the  same  place.  Subroutine  Mwfm  grid  in  Appendix  E 
provides  the  FORTRAN  code  used  to  compute  the  grid  point’s  mountain  wave  values. 

3.4.2  PIREP  Index  Assignment  Implementation.  To  calibrate  the  index  values  to 
turbulence  intensity,  pilot  report  data  and  the  model  analyses  were  used.  The  PIREPS 
were  provided  by  AECCC  for  the  period  of  May  20,  2003  to  October  15,  2003.  A  large 
file  containing  all  the  PIREPS  from  this  period  was  parsed  by  a  shell  script  to  only 
include  the  PIREPS  in  a  six  hour  window  centered  on  the  valid  time  of  the  forecast. 

Then,  for  each  pilot  report,  the  levels  above  and  below  the  PIREPS  were  determined 
based  on  the  level  height  for  the  Ellrod-2  index.  Einally,  the  maximum  of  each  index  of 
the  eight  grid  points  surrounding  the  location  of  the  pilot  report  were  assigned  to  each 
PIREP.  This  method  matched  the  method  used  in  the  Brown  et  al.  (2000)  study.  Refer  to 
Appendix  E  subroutine  PR  int  for  the  EORTRAN  code  used  to  assign  index  values  to 
PIREPS. 
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3.4.3  Composite  Index  Implementation.  The  composite  index  was  computed  for  all  the 
levels  and  layers  in  the  study.  To  determine  it,  the  FORTRAN  program  computed  a 
value  of  0,  1,  2,  or  3  for  none,  light,  moderate,  or  severe  turbulence  for  each  grid  point  in 
the  array.  Turbulence  values  were  assigned  by  comparing  each  index  to  each  threshold, 
as  described  in  Section  4.2,  for  every  layer  or  level.  When  three  adjacent  layers  or  levels 
exceeded  an  index’s  threshold,  and  that  grid  point  had  not  already  been  assigned  a  greater 
turbulence  value,  that  grid  point  was  assigned  the  corresponding  turbulence  value.  For 
example,  if  the  450  mb,  500  mb,  and  550  mb  Ellrod-2  index  values  all  exceeded  the 
moderate  threshold  for  a  grid  point  and  the  grid  point  had  previously  been  assigned  a  0 
(none)  or  1  (light),  the  grid  point  was  reassigned  a  2  (moderate).  If  the  same  grid  point 
had  already  been  assigned  a  3  (severe),  it  would  have  remained  a  3.  Refer  to  Appendix 
E,  subroutine  Comp  make  for  the  EORTRAN  code  used  to  create  this  chart. 

Additionally,  the  maximum  value  for  each  index  for  each  grid  point  was  also  computed 
and  stored  in  an  array.  Einally,  these  four  additional  fields  (composite,  max  Ellrod-2, 
max  MWEM  VI  .1,  and  max  MWEM  V2.1)  were  stored  at  the  end  of  the  data  file  with 
the  individual  indices  for  display  using  GrADS  software. 

3.4.4  GrADS  Implementation.  GrADS  software  was  used  to  make  the  necessary  image 
files  for  the  forecaster  graphics.  GrADS  works  with  unformatted  EORTRAN  data  stored 
in  four  dimensions  very  easily  (longitude,  latitude,  height,  time).  The  data  fdes  created 
by  this  EORTRAN  program  saved  the  index  values  in  three  dimensions  (longitude, 
latitude,  height).  The  longitude  array  elements  were  stored  from  140W-30W  and  the 
latitude  array  elements  were  stored  from  20N-60N.  The  height  elements  were  stored 
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such  that  the  layers  from  100-150  mb  to  700-750  mb  of  the  MWFM  indices  corresponded 
to  the  levels  from  150  mb  to  750  mb  for  the  Ellrod-2  index.  The  composite  layer  was 
stored  as  an  additional  height  layer  after  the  bottom  layer  in  the  array.  GrADS  relies  on 
metadata  files  to  tell  it  how  to  interpret  the  data  it  reads.  GrADS  also  has  a  scripting 
feature  to  string  together  a  sequence  of  commands.  The  metadata  file  and  script  used  to 
generate  the  imagery  are  provided  in  Appendix  C.  Information  on  using  GrADS  was 
obtained  from  on-line  documentation  (GrADS  2003). 


lnitial_Dote=Q31 106_0Q,  Fcst_tim€=  I  2Z,  COMPOSITE  ANALYSIS 


Figure  5.  Example  composite  chart 


Eor  viewing  the  data,  the  indices  were  contoured  on  a  map  of  the  CONTIS.  The 
contours  were  determined  by  the  thresholds  as  determined  by  the  PIREP  data.  The 
threshold  values,  as  determined  in  Section  4.2,  were  2.16,  3.00,  and  8.37  for  the  Ellrod-2 
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index,  0.00033,  0.00140,  and  undefined  for  MWFM  V 1 . 1  index,  and  0.012,  0.020,  and 
0.036  for  the  MWFM  V2.1  index.  The  eharts  developed  by  GrADS  were  saved  in  “Joint 
Photographie  Experts  Group”  (JPEG)  format. 

To  reduee  the  work  required  by  the  foreeaster,  a  eomposite  image  was  ereated 
from  the  eomposite  index.  This  image  shaded  eaeh  respeetive  threshold  with  values  of  2 
and  3  being  depieted  as  gray  shades  in  figure  5.  In  figure  5,  regions  of  moderate 
turbulenee  are  shaded  lighter  than  regions  of  severe  turbulenee  and  are  deseribed  in 
Seetion  3.4.3.  The  intent  was  to  allow  the  foreeaster  to  traee  the  regions  outlined  by  the 
eomposite  ehart  and  use  the  layer  viewer  to  determine  the  levels  of  the  atmosphere  that 
indieate  the  potential  for  CAT.  This  image  was  saved  in  JPEG  format  also. 

3.4.5  Turbulence  Layer  Viewer  Implementation.  The  turbulenee  layer  viewer  provided  a 
quiek  method  of  interrogating  the  atmosphere  and  determining  the  turbulent  levels. 
Beeause  of  the  relatively  large  size  an  image  must  be  to  properly  interrogate  it  with  a 
eomputer,  foreeasters  are  ordinarily  limited  to  examining  eaeh  layer  individually.  The 
layer  viewer  eireumvents  that  problem  by  allowing  foreeasters  to  view  the  same  portion 
of  several  layers  of  the  atmosphere  at  the  same  time. 

The  layer  viewer  requires  frames  and  javaseript  to  work  and  was  designed  to 
work  with  the  Mierosoft®  Internet  Explorer  V6.0  web  browser.  It  is  eomposed  of  two 
primary  areas.  The  first  area  holds  the  layers.  To  speed  up  the  transition  between  layers, 
the  viewer  immediately  loads  all  the  available  layers  into  memory  despite  only  displaying 
a  few  of  them.  In  this  study,  four  layers  of  the  14  were  displayed  in  the  first  area  but  that 
number  eould  be  easily  modified.  The  seeond  primary  area  holds  the  two  tools  needed 
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to  control  the  layers.  The  first  tool  is  the  up/down  buttons.  Clicking  on  either  of  these 
buttons  shifts  the  four  displayed  layers  up  or  down  by  one  layer.  The  other  tool  is  a 
vertical  scrollbar.  Because  of  the  viewing  size  of  the  four  images,  only  a  portion  of  each 
image  in  the  first  area  display  region  is  viewable  at  a  time.  Adjusting  the  scrollbar  shifts 
the  viewable  region  in  each  window  of  the  first  area  display  region  to  the  same  place. 
This  feature  may  not  work  properly  with  a  Netscape®  browser  or  older  Internet  Explorer 
versions.  The  images  used  in  this  study  did  not  require  a  horizontal  scrollbar.  If  a 
horizontal  scrollbar  were  required,  a  few  simple  modifications  of  the  javascript  used  by 
the  vertical  scrollbar  are  all  that  would  be  required.  Refer  to  Figure  6  for  an  example  of 
the  layer  viewer  and  Appendix  D  for  a  summary  of  the  html  and  javascript  code  used  in 
making  the  layer  viewer. 
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Figure  6.  Turbulence  layer  viewer 
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Aside  from  the  two  primary  areas,  additional  areas  provide  an  approximate  flight 
level  for  each  layer  and  a  title.  The  approximate  flight  level  was  calculated  for  layers 
below  the  tropopause  layer  using  equations  14  and  15. 
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Equation  14  is  based  on  the  U.S.  Standard  Atmosphere  below  the  approximate 
tropopause  (8  km  or  about  the  200-250  mb  layer).  In  this  equation,  z  is  the  height,  is 

the  dry  air  gas  constant  287.04  J ■kg~^  ■  K  \  T  is  the  average  temperature  288.15  K,  g  is 
the  acceleration  due  to  gravity  9.80665  m-s~^ ,  is  the  reference  pressure  1000  mb,  and 


p  is  the  pressure  at  the  level  of  interest.  Equation  15  averages  the  height  for  the  top  and 


bottom  level  of  each  layer. 


Table  5,  List 

:  of  flight  levels  based  on  pressure 

Layer 

Flight  Level 

Layer 

Flight  Level 

100-150  mb 

FL480 

450-500  mb 

FL200 

150-200  mb 

FL420 

500-550  mb 

FL170 

200-250  mb 

FL360 

550-600  mb 

FL150 

250-300  mb 

FL320 

600-650  mb 

FL130 

300-350  mb 

FL280 

650-700  mb 

FL110 

350-400  mb 

FL250 

700-750  mb 

FL090 

400-450  mb 

FL220 

Eor  each  of  the  layers,  the  height  was  verified  with  a  sample  of  the  model  output. 
Additional  comparison  was  performed  visually  with  a  skew-T  of  the  El.S.  standard 
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atmosphere.  It  is  important  to  remember  these  heights  are  estimates  of  the  heights  in  the 
atmosphere  and  are  to  be  used  as  a  guide.  The  values  determined  using  these  methods 
are  presented  in  Table  5. 

The  region  of  the  atmosphere  ehosen  for  the  upper  level  turbulenee  eharts  poses  a 
eoneern  to  this  method  of  height  determination.  When  flying  at  or  above  FL180,  aviators 
set  their  altimeter  to  the  standard  surfaee  pressure,  29.92  inehes  Hg.  When  flying  below 
FL180,  aviators  use  a  nearby  station  pressure  (FAA  2003).  This  differenee  in  altimeter 
settings  ean  produee  differenees  in  indieated  altitude  of  several  hundred  feet.  For  the 
purpose  of  this  study,  the  potential  differenee  in  altitude  was  eonsidered  as  a  possible 
souree  for  error  but  no  efforts  were  made  to  aeeount  for  it  sinee  the  eharts  assign  a  region 
of  altitudes  as  potentially  turbulent  levels  and  not  speeifie  heights. 

The  data  file  eontaining  the  three  index  values  and  eomposite  data  was  used  to 
populate  the  turbulence  layer  viewer.  The  images  were  generated  using  GrADS  scripts  as 
described  in  section  3.4.4. 

3.5  Turbulence  Chart  Methodology 

After  determining  the  thresholds,  the  test  method  was  verified  using  PIREPS. 
Producing  the  turbulence  charts  with  the  test  method  required  two  steps.  First,  depending 
on  the  desired  minimum  turbulence,  the  contour  indicating  the  desired  turbulence 
intensity  was  outlined.  For  this  study,  MOG  turbulence  intensity  was  chosen  since  as  a 
matter  of  practice,  weather  squadrons  only  forecast  regions  of  moderate  or  greater 
turbulence.  When  two  or  more  areas  suggesting  MOG  turbulence  were  very  close 
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together,  one  eirele  eneompassing  all  the  areas  was  drawn.  Seeond,  using  the  layer 
viewer,  the  altitudes  where  turbulenee  was  present  were  loeated  and  the  eireled  regions 
on  the  map  were  labeled  with  the  estimated  base  and  top  of  suspeeted  turbulenee.  This 
oeeasionally  required  dividing  large  outlined  regions  into  parts  to  make  a  more 
representative  visualization. 

Despite  every  attempt  to  make  a  purely  objeetive  foreeast  tool,  several  aspeets  of 
this  method  required  foreeaster  subjeetivity.  First,  when  two  or  more  areas  on  the 
eomposite  ehart  suggesting  MOG  turbulenee  were  very  elose  together,  one  eirele 
eneompassing  all  the  areas  was  drawn.  This  improved  the  appearanee  of  the  charts  but 
required  the  forecaster  to  determine  when  two  regions  were  close  enough  to  be  joined. 
Second,  there  were  several  situations  that  often  arose  when  using  the  layer  viewer  which 
required  some  interpretation.  First,  within  a  layer,  different  intensities  were  present 
within  one  of  the  highlighted  regions.  In  general,  the  level  was  included  if  more  than  half 
of  the  layer  contained  any  indication  of  turbulence  with  at  least  some  of  that  meeting  the 
moderate  or  greater  threshold.  The  other  situation  involved  interpreting  between  two  or 
more  levels.  Occasionally,  two  levels  which  obviously  needed  to  be  included  were 
divided  by  a  layer  that  appeared  less  likely  to  contain  turbulence.  In  this  instance,  any 
indication  of  moderate  or  greater  turbulence  by  the  intermediate  layer  led  to  inclusion  of 
that  altitude.  Furthermore,  the  top  or  bottom  of  the  turbulent  region  rarely  stopped  on  a 
layer  exactly.  Determining  the  actual  heights  of  the  top  and  bottom  required  some 
interpretation  by  the  forecaster.  All  of  the  subjective  decisions  were  made  as  a  matter  of 
practicality  and  are  similar  in  practice  to  steps  used  by  USAF  weather  forecasters  to 


41 


estimate  heights  and  regions  of  turbulenee.  These  simplifications  allow  the  aviators  and 
mission  planners  to  review  a  usable,  easily  readable  turbulence  chart. 

Verification  of  charts  produced  by  the  test  method  and  the  weather  squadrons 
were  also  done  using  the  same  PIREPS  database.  The  results  of  both  methods  using  the 
statistics  discussed  in  Section  3.6.2  are  detailed  in  Chapter  4. 

3.6  Statistical  Methods 

The  statistical  analysis  performed  during  this  research  can  be  divided  into  two 
primary  parts.  First,  the  thresholds  for  each  of  the  three  indices  needed  to  be  determined. 
This  was  done  to  provide  a  basis  for  the  forecast  tools.  The  primary  goal  of  this  stage 
was  to  determine  usable  values  of  each  index  that  correspond  with  light,  moderate,  and 
severe  turbulence.  These  thresholds  were  then  used  to  create  the  images  used  by  the 
turbulence  viewer  software  to  assist  the  forecaster  in  turbulence  chart  development.  The 
second  stage  of  the  statistical  analysis  involved  verifying  forecasts  produced  using  the 
test  method  with  PIREPS  and  comparing  it  to  verification  of  actual  forecasts  produced  by 
USAF  weather  squadrons. 

3.6.1  Threshold  Determination.  The  thresholds  were  determined  with  a  FORTRAN 
program  which  processed  every  possible  combination  of  indices  looking  for  the  ideal 
combination  based  on  contingency  table  statistics  as  described  in  Table  1.  Several 
different  criterion  were  compared  to  determine  the  thresholds  including  the  PODy 
maximum,  PODn  maximum,  FAR  minimum,  CSI  maximum,  PC  maximum,  TSS 
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maximum,  closest  BIAS  score  to  one,  the  first  instance  when  PODy  was  less  than  or 
equal  to  PODn,  the  maximum  HSS,  and  the  maximum  GSS. 

The  HSS  and  GSS  are  statistics  developed  to  determine  the  actual  skill  of  a 
forecast  relative  to  random  forecasts.  According  to  Wilks  (1995),  the  HSS  is  the  most 
common  skill  score  used  to  summarize  square  contingency  tables.  It  essentially 
compares  the  forecasts  to  a  generic  reference  forecast,  using  the  equation  in  Table  1,  with 
the  same  marginal  distributions  of  forecasts  and  observations  as  the  actual  contingency 
table.  The  HSS  returns  a  value  less  than  or  equal  to  one.  A  score  of  one  signifies  a 
perfect  forecast  while  forecasts  worse  than  the  generie  referenee  foreeast  receive  scores 
less  than  zero.  Schaefer  (1990),  describes  the  GSS  as  the  CSI  corrected  for  the  number  of 
hits  expected  by  chance.  A  perfect  forecast  receives  a  score  of  one.  The  closer  a 
forecast  is  to  being  purely  random,  the  closer  the  GSS  will  be  to  zero  (Stephenson  2000). 
The  equations  used  in  caleulating  the  HSS  and  GSS  are  listed  in  Table  1. 

The  FORTRAN  program  contained  four  nested  loops.  The  loops  cycled  through 
every  combination  of  the  three  indices  over  user-determined  ranges  of  values  at  user- 
determined  inerements.  Additionally,  the  program  could  be  configured  to  analyze 
different  flight  levels,  turbulence  intensities,  latitude/longitude  ranges,  and  different 
lengths  of  time  centered  on  the  foreeast  hour  for  the  same  statistics.  After  eomputing  the 
statistics,  the  FORTRAN  program  saved  the  results  to  an  output  file  for  review. 

After  determining  the  thresholds,  the  jl  t^st,  as  described  in  Montgomery  and 
Runger  (2003),  was  performed  on  the  contingency  tables  to  test  for  independence.  In  this 
test,  the  null  hypothesis  was  that  forecasts  were  independent  of  observations.  An  alpha  of 
0.05  was  used  to  reject  the  null  hypothesis  and  determine  whether  forecasts  and 
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observations  of  turbulence  were  dependent  on  each  other.  The  test  statistic  in  a  test  is 
computed  using  the  equation: 


I  J 


(16) 


where  O^.  is  the  observed  contingency  table  values,  E^j  is  the  expected  table  values,  and  i 

and  j  are  the  rows  and  columns.  The  observed  contingency  table  is  the  contingency  table 
calculated  from  the  data  while  the  expected  contingency  table  is  computed  using: 


tot  .{tot  j) 
total 
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where  tot.  represents  the  total  of  row  i  in  the  observed  contingency  table,  tot.  is  the  total 


of  column  j  in  the  observed  contingency  table,  and  total  is  the  sum  of  all  the  cells  in  the 
observed  contingency  table.  After  determining  the  yl  value  for  the  data,  it  was  compared 
to  the  table  value  of  yl  for  an  alpha  of  0.05  with  1  degree  of  freedom  or  3.84.  If  the 
computed  yl  value  exceeded  3.84,  the  null  hypothesis  was  rejected. 


3.6.2  PIREP  Verification.  After  the  thresholds  were  determined,  the  layer  viewer 
software  was  configured  with  the  thresholds  to  allow  for  forecasts  to  be  developed. 

These  forecasts  were  made  over  the  CONUS  from  October  21,  2003  to  October  27,  2003 
and  from  November  5,  2003  to  November  11,  2003.  Contingency  tables  combining 
PIREPS  and  forecasts  were  used  to  compute  the  same  statistics  discussed  in  section  3.6.1. 
Additionally,  forecasts  produced  by  USAF  weather  squadrons  over  the  CONUS  for  the 
same  dates  were  verified  using  the  same  methods.  The  USAF  weather  squadron  forecasts 
originated  from  stitched  charts  produced  by  the  26**’  OWS  at  Barksdale  AFB.  These 
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stitched  charts  represented  the  combined  forecasts  of  turbulence  charts  produced  by  the 


15*  OWS  at  Scott  AFB,  25*  OWS  at  Davis-Monthan  AFB,  26*  OWS  at  Barksdale  AFB, 
and  28*  OWS  at  Shaw  AFB. 
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IV.  Results 


4.1.  Introduction 

This  chapter  summarizes  the  results  from  the  various  eontingeney  table  analyses 
performed  using  the  PIREPS  and  the  various  foreeasts.  It  begins  with  a  diseussion  of  the 
various  analyses  performed  on  the  data  to  determine  the  different  thresholds.  This 
ineludes  testing  for  whether  the  thresholds  vary  by  altitude  or  geographie  region.  It  also 
ineludes  testing  the  thresholds  for  turbulenee  greater  than  or  equal  to  light,  greater  than  or 
equal  to  moderate,  and  greater  than  or  equal  to  severe.  After  eompleting  the  threshold 
determination,  the  established  thresholds  were  entered  into  the  FORTRAN  program  that 
eomputed  the  indiees  and  the  program  was  rerun  to  generate  the  produets  used  by  the  test 
foreeast  method.  The  seeond  seetion  deseribes  the  results  of  the  foreeast  eomparison 
between  the  test  method  and  operational  foreeasts.  It  foeuses  on  the  relative  aeeuraey  of 
the  test  method  to  the  foreeasts  eurrently  produeed  by  the  various  weather  squadrons. 

4.2.  Threshold  Determination 

Thresholds  were  determined  using  the  subset  of  PIREPS  within  six  hours  of  the 
model  analysis  runs  from  May  20,  2003  to  October  15,2003.  For  example,  for  a  model 
foreeast  valid  at  12Z,  any  PIREP  between  0900  and  1459  were  used.  Elsing  this  subset  of 
PIREPS  provided  a  sample  size  of  15,745.  In  determining  the  thresholds,  the  Ellrod-2 
index  threshold  was  found  for  eaeh  intensity  level  first.  This  was  done  by  setting  the 
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ranges  for  the  MWFM  indiees  with  a  value  of  99.0.  This  value  was  well  outside  of 
maximum  value  of  either  MWFM  index.  After  determining  this  threshold,  the  MWFM 
V2.1  threshold  was  found  with  the  Ellrod-2  index  set  to  its  already  determined  threshold. 
Finally,  the  MWFM  Vl.l  threshold  was  found  with  the  other  two  indiees  set  to  their 
established  thresholds.  This  order  of  threshold  determination  was  ehosen  beeause  the 


MWFM  Vl.l  threshold  tends  to  eontain  output  for  the  region  direetly  above  a  ridge  while 
MWFM  V2.1  propagates  energy  downstream  making  it  provide  output  over  larger 
regions  of  the  map.  The  order  intended  to  ensure  the  best  ehanee  that  eaeh  index 
eontributed  to  finding  positive  reports  of  turbulenee  the  other  indiees,  eovering  larger 
areas,  missed.  The  eriteria  most  often  used  to  set  the  thresholds  ineluded  the  skill  seores 


and  the  point  where  PODy  erossed  PODn.  The  results  of  threshold  determination  are 


provided  in  Table  6. 


Table  6,  Thresholds  determined  for  indices 


Light 

Moderate 

Severe  | 

Chart 

Index 

TH 

PODy 

PODn 

TH 

PODy 

PODn 

TH 

PODy 

PODn 

Layers 

Ellrod  2 

2.16 

0.50 

0.63 

3.00 

0.45 

0.67 

8.37 

0.21 

0.94 

MWFM  V1 

0.00033 

0.67 

0.49 

0.0014 

0.56 

0.60 

- 

- 

- 

MWFM  V2 

0.012 

0.66 

0.49 

0.02000 

0.54 

0.61 

0.03600 

0.33 

0.87 

From  Table  6,  the  light  and  moderate  thresholds  were  quite  elose  to  eaeh  other  for 
eaeh  of  the  indiees.  This  would  make  separating  regions  of  light  and  moderate 
turbulenee  diffieult.  Sinee  Air  Foree  weather  squadrons  only  prediet  regions  of  moderate 
or  greater  turbulenee  in  praetiee,  this  problem  should  not  have  signifieantly  affeeted  the 
skill  of  the  test  foreeast  method. 

Threshold  determination  also  provided  several  noteworthy  results  not  listed  in 
Table  6.  First,  there  were  only  24  instanees  of  severe  turbulenee  in  this  sample  making 
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determination  of  the  thresholds  for  that  intensity  unreliable.  This  led  to  the  establishment 
of  an  almost  arbitrary  value  for  the  severe  turbulenee  thresholds  and  no  value  being 
determined  for  the  MWFM  Vl.l  index  severe  threshold.  As  sueh,  any  indieation  of 
severe  turbulenee  by  the  foreeast  method  was  ignored  for  this  study  and  only  moderate 
regions  were  highlighted  for  the  verifieation  study  of  Seetion  4.3.  Additionally,  the  'i 
test  was  performed  on  the  eontingeney  table  produced  by  the  combined  thresholds  for 
moderate  turbulence.  The  resulting  value  of  193.5323  suggested  there  was  a  5.39x10'  % 
chance  of  incorrectly  rejecting  the  null  hypothesis  with  this  y  value.  This  provided 
strong  evidence  to  reject  the  null  hypothesis  and  conclude  the  index  values  offer  some 
skill  in  predicting  turbulence.  Another  interesting  fact  of  this  test  was  that  FAR  did  not 
increase  with  the  addition  of  MWFM  indices  to  the  forecast.  In  fact  FAR  decreased  from 
0.8394  to  0.8339  to  0.8338  as  the  MWFM  indices  were  added  to  Ellrod-2  index  when 
computing  the  moderate  layer  threshold.  This  provides  another  indication  that  the 
MWFM  indices  improve  the  forecast.  Additionally,  the  MWFM  indices  provide  a  larger 
improvement  in  the  PODy  than  the  degradation  in  the  PODn,  suggesting  they  add  value 
to  the  forecast. 

After  determining  the  thresholds,  it  was  necessary  to  see  if  they  change 
significantly  for  different  geographic  regions  or  altitudes.  This  test  was  performed  by 
setting  the  indices  to  the  established  thresholds  and  calculating  the  statistics  over  different 
regions  or  altitudes.  For  this  test,  there  were  nine  different  regions  tested  for  moderate 
turbulence.  The  regions  included  the  following  areas:  A)  Entire  region,  B)  EEIOO- 
EE200,  C)  EE200-EE300,  D)  EE300-EE400,  E)  EE400-EE500,  E)  1 10W-70W  and  20N- 
40N,  G)  1 10W-70W  and  40N-60N,  H)  70W-30W  and  20N-40N,  and  I)  70W-30W  and 


48 


40N-60N.  The  resulting  statistics  are  listed  in  Table  7.  The  statistics  varied  considerably 
for  different  regions  and  levels  in  the  atmosphere  suggesting  different  thresholds  based 
on  geographic  region  or  altitude  might  produce  more  accurate  results.  In  addition, 
several  of  these  samples  failed  the  x  test  for  independence.  Those  tests  include  E,  H,  and 
I  with  'i  values  of  2.00,  1 .38,  and  3.44  respectively.  Each  of  these  values,  with  an  a  of 
0.05,  fail  to  reject  the  null  hypothesis.  In  test  H,  the  “yes/yes”  cell  only  expected  1.62 
occurrences.  According  to  Wilks  (1995),  cells  containing  less  than  five  instances  of  an 
expected  event  should  be  avoided.  Also,  in  tests  H  and  I,  there  were  only  13  and  35 
positive  reports  of  moderate  or  greater  turbulence.  This  lack  of  positive  reports  decreased 
the  credibility  of  the  statistics  listed  for  those  rows.  Despite  the  lack  of  validity  regarding 
these  two  regions,  there  does  appear  to  be  some  regional  variation  in  the  thresholds. 


Table  7.  Statistics  for  different  regions 


Test 

PODy 

PODn 

FAR 

CSI 

PC 

BIAS 

TSS 

HSS 

GSS 

A 

15745 

0.5628 

0.5983 

0.8225 

0.1560 

0.5935 

3.1703 

0.1611 

0.0840 

0.0438 

B 

2136 

0.5367 

0.6001 

0.6004 

0.2971 

0.5791 

1 .3432 

0.1369 

0.1260 

0.0672 

C 

1161 

0.6187 

0.5102 

0.6240 

0.3053 

0.5452 

1 .6453 

0.1288 

0.1102 

0.0583 

D 

11987 

0.5696 

0.6030 

0.8911 

0.1006 

0.6004 

5.2295 

0.1726 

0.0588 

0.0303 

E 

463 

0.4487 

0.6364 

0.8000 

0.1606 

0.6048 

2.2436 

0.0851 

0.0569 

0.0293 

F 

4837 

0.5555 

0.6466 

0.6498 

0.2735 

0.6233 

1.5862 

0.2021 

0.1694 

0.0925 

G 

4922 

0.5822 

0.5503 

0.7913 

0.1815 

0.5557 

2.7899 

0.1325 

0.0774 

0.0402 

H 

869 

0.2308 

0.8773 

0.9722 

0.0254 

0.8677 

8.3077 

0.1081 

0.0235 

0.0119 

1 

5167 

0.6000 

0.5563 

0.9909 

0.0091 

0.5566 

65.6571 

0.1563 

0.0047 

0.0024 

4.S  Forecast  Comparison 


To  compare  forecast  skill,  12-hour  and  24-hour  forecasts  were  made  using  the  test 


method  and  verified  with  PIREPS.  Additionally,  12-hour  and  24-hour  operational 


forecasts  produced  by  the  Air  Eorce  were  also  verified  with  PIREPS.  When  initially 
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attempting  to  implement  the  test  method  as  an  operational  forecast  tool,  the  indices 
produced  output  that  was  very  confusing  and  that  would  have  been  difficult  to  use 
operationally.  The  composite  images  produced  with  these  thresholds  were  very 
congested  with  contours,  as  depicted  by  the  image  on  the  left  in  Figure  7.  As  a  result,  the 
composite  image  only  depicted  regions  with  turbulence  when  three  adjacent  layers 
indicated  turbulence  at  a  grid  point  as  described  in  section  3.4.3.  This  modification 
produced  the  image  on  the  right  in  Figure  7  which  was  much  easier  to  contour. 


lime-l2Z,  COMPOSITE  ANALYSIS 


Ar'^  ^ 

s' 


V«tto««modero1e  (tfituience  Red*S«vere  lurbui«nc«  {pos»ble) 

Figure  7.  Examples  of  composite  level  output 


After  adjusting  the  method  for  computing  the  composite  chart,  forecasts  were 
made  using  the  test  method.  Typically,  these  hand-plotted  forecasts  took  less  than  five 
minutes  to  produce.  Next,  these  test  forecasts  and  actual  forecasts  produced  by 
operational  weather  squadrons  were  verified  for  accuracy  with  PIREPS  from  October  21, 
2003  to  October  27,  2003  and  from  November  5,  2003  to  November  12,  2003.  To  verify, 
PIREPS  were  compared  by  location  and  height  to  the  forecasts.  In  order  to  be  counted, 
PIREPS  needed  to  be  above  the  CONUS  since  the  weather  squadrons  did  not  forecast 
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outside  the  CONUS.  When  a  PIREP  fell  within  a  region  highlighted  to  suggest 
turbulenee,  it  was  eonsidered  foreeasted.  Only  PIREPS  of  moderate  or  greater  intensity 
were  eonsidered  observed. 

The  PIREPS  were  plotted  using  GrADS  on  a  ehart  that  eontained  the  test  method 

foreeast  regions  so  the  eomparison  was  quite  accurate  for  the  test  method.  Operational 

forecasts  were  compared  with  PIREPS  using  estimation.  To  ensure  the  weather  squadron 

forecast  skill  was  maximized,  PIREPS  appearing  very  close  to  a  region  of  forecasted 

turbulence  were  given  the  benefit  of  the  doubt.  In  general,  these  events  were  scored  as 

either  “yes/yes”  or  “no/no”.  Furthermore,  this  study  only  predicted  regions  of  moderate 

or  greater  turbulence.  There  were  many  PIREPS  of  light  turbulence  reported  in 

forecasted  regions  of  moderate  turbulence.  Those  regions  were  scored  as  forecasted 

“yes”  but  observed  “no”  or  “yes/no”.  The  test  method  and  operational  forecast 

contingency  tables  are  provided  in  tables  8  and  9.  The  x  test  scores  of  189.9  for  the  test 

method  and  35.97  for  the  operational  forecasts  are  sufficient  to  reject  the  null  hypothesis 

leading  to  the  conclusion  that  forecasts  and  observations  are  not  independent. 

Table  8.  Test  method  contingency  table  Table  9.  OWS  forecast  contingency  table 
Observed  Observed 


Yes 

No 

Yes 

No 

Forecasted 

Forecasted 

Yes 

325 

478 

Yes 

87 

130 

No 

246 

1371 

No 

484 

1719 
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These  tables  suggest  signifieant  differenees  exist  in  the  two  methods.  In  general, 
the  operational  foreeasts  tend  to  significantly  under  forecast  turbulence  event  relative  to 
the  test  method,  as  evidenced  by  the  large  number  in  the  “no/yes”  cell.  At  the  same  time, 
the  test  method  appears  to  significantly  over  forecast  turbulence.  This  apparent  weakness 
is  tempered  by  the  short-lived  nature  of  CAT  and  the  low  reliability  of  the  “yes/no”  cell 
as  a  determinant  of  forecast  skill.  Also,  negative  reports  of  turbulence  would  often 
appear  very  close  to  positive  reports  of  turbulence  in  time  and  location  further  validating 
the  claim  that  CAT  is  short-lived  and  localized.  Also,  the  operational  forecasts  did  not 
appear  to  consider  turbulence  caused  by  mountain  waves.  There  was  not  one  region 
depicted  on  the  charts  over  these  two  weeks  where  the  weather  squadrons  used  solid  blue 
lines  to  depict  mountain  wave  turbulence.  This  may  be  due  to  mislabeling  since  there 
were  several  regions  which  suggested  the  forecaster  considered  the  turbulence  to  be 
mountain  wave  turbulence  although  they  used  the  general  CAT  line  style  to  circle  them. 


Table  10.  Verificai 

tion  $tai 

tistics 

PODy 

PODn 

FAR 

CSI 

PC 

BIAS 

TSS 

HSS 

GSS 

Test  Method 

0.569 

0.741 

0.595 

0.310 

0.701 

1.406 

0.311 

0.272 

0.231 

Operational 

Forecasts 

0.152 

0.930 

0.599 

0.124 

0.746 

0.380 

0.082 

0.104 

0.114 

Instances  where  the  test  method  did  not  forecast  regions  that  observed  turbulence 
were  explored  further  to  determine  the  cause.  There  were  several  possible  reasons  the 
test  method  did  not  highlight  these  regions  of  turbulence.  Several  observed  instances  of 
turbulence  appear  to  have  been  caused  by  convection.  These  instances  would  not  meet 
the  definition  of  CAT  but  were  left  in  the  tables  since  this  is  a  comparative  study  of 
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forecast  methods  and,  in  general,  neither  method  forecasted  those  regions.  The  most 
common  cause  of  missed  regions  of  turbulence  was  near  misses  where  the  turbulence 
appeared  just  outside  of  the  highlighted  areas.  Many  of  these  near  misses  were  captured 
with  the  actual  layer  viewer.  Table  10  provides  the  detailed  breakdown  of  the  statistics 
produced  from  the  contingency  tables  for  the  test  method  and  operational  forecasts. 

The  test  method  appears  superior  to  operational  forecasts.  The  test  method  scored 
significantly  higher  PODy,  TSS,  HSS,  and  GSS  over  the  operational  forecasts.  The 
higher  PODn  of  the  operational  forecasts  is  due  in  large  part  to  the  significantly  lower 
volume  coverage  of  the  forecasts.  The  bias  scores  suggest  the  test  method  tends  to  over¬ 
forecast  turbulence  while  the  operational  forecast  tends  to  under-forecast  turbulence. 
However,  the  test  method  bias  is  closer  to  one  suggesting  it  tends  to  over-forecast  less 
than  the  operational  forecast  tends  to  under-forecast. 
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V.  Conclusions  and  Recommendations 


5.1  Conclusions 

For  many  years,  meteorologists  have  attempted  to  prediet  oeeurrenees  of  elear-air 
turbulenee.  In  reeent  years,  various  numerieal  indiees  have  been  developed  and  tested  to 
quantitatively  prediet  CAT.  It  had  been  shown  that  numerieal  models  outperform 
eoneeptual  models  when  foreeasting  CAT  (Dutton  1980).  This  study  seems  to  baek  up 
that  assertion.  This  researeh  foeused  on  developing  a  new  teehnique  for  predieting  CAT 
that  ineorporated  three  numerieal  model  indiees,  the  Ellrod-2  index  and  the  two  MWFM 
indiees.  Ineluding  tools  for  predieting  both  primary  eauses  of  CAT  provided  a 
eomprehensive  predietion  method.  The  teehnique  was  developed  to  meet  the  needs  of 
Air  Foree  operational  weather  squadron  foreeasters  for  a  eomprehensive,  easy-to-use 
automated  foreeast  method.  Furthermore,  the  teehnique  is  adaptable  to  indiees  other  than 
those  ehosen  in  this  study. 

The  mountain  wave  foreeast  model  indiees  showed  promise  in  predicting  CAT 
regions  in  the  troposphere.  The  V2.1  index  provides  more  comprehensive  coverage  than 
the  Vl.l  index.  However,  MWFM  V2.1  does  not  appear  to  capture  every  region  of 
mountain  wave  turbulence,  and  the  Vl.l  index  appears  to  add  some  value  to  the  overall 
forecast.  Although  most  turbulence  encounters  appear  to  be  represented  by  an  index 
depicting  Kelvin-Helmholtz  instabilities,  there  were  several  instances  that  the  MWFM 
captured  but  that  the  Ellrod-2  index  missed.  In  this  study,  the  mountain  wave  indices 
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appeared  to  capture  about  one  in  every  ten  correctly  forecasted  moderate  or  greater 
turbulence  encounters  throughout  the  column,  not  just  near  the  tropopause. 

The  test  method  appears  to  provide  superior  forecasts  to  those  currently  produced 
by  the  operational  weather  squadrons  who  appear  to  primarily  use  conceptual  models  in 
their  turbulence  chart  development.  This  is  most  likely  the  cause  of  the  significant  under¬ 
forecasting  and  inaccuracy  in  the  weather  squadron  products  used  for  this  research. 

While  the  test  method  provided  a  large  number  of  near  misses  on  occurrences  of 
turbulence,  PIREPS  of  turbulence  were  often  quite  far  from  the  nearest  turbulent  region 
indicated  by  the  weather  squadron  forecasts.  However,  the  test  method  charts  tend  to 
encompass  a  large  portion  of  the  atmosphere  making  turbulent  conditions  appear  quite 
widespread.  Furthermore,  the  test  method  appears  to  miss  quite  a  few  instances  of 
turbulence  despite  the  layer  viewer  depicting  those  regions  as  turbulent.  These  misses  are 
caused  by  the  implementation  of  the  composite  chart  as  described  in  Section  3.4.3.  Also, 
because  of  the  significant  coverage  of  turbulent  regions  by  the  test  method,  using  the  test 
method  is  not  as  straightforward  as  initially  thought.  This  was  especially  true  over  the 
Rocky  Mountains  as  they  often  have  numerous  levels  and  layers  where  the  mountain 
wave  forecast  models  predicted  turbulence. 

This  study  only  focused  on  moderate  or  greater  turbulence.  The  length  of  time 
used  in  this  study  did  not  provide  an  adequate  number  of  severe  PIREPS  (due  to  their 
limited  availability)  to  properly  threshold  severe  turbulence  events.  The  few  instances  of 
severe  turbulence  actually  depicted  in  this  study  did  allow  for  some  degree  of  threshold 
determination;  however,  the  results  are  unreliable  and,  therefore,  deemed  inconclusive. 
Properly  thresholding  severe  turbulence  would  require  a  much  longer  test  period. 
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encompassing  many  more  severe  turbulence  encounters  than  identified  in  the  sample 
analysis. 

Overall,  the  test  method  appears  to  offer  promise  as  a  tool  that  simplifies  an 
otherwise  complex  process.  The  time  required  to  develop  forecasts  for  the  CONUS  using 
this  new  tool  should  be  less  than  five  minutes  per  forecast  hour  shortly  after  the 
forecaster  masters  the  technique  and  the  OPS  II  drawing  features.  Furthermore,  the 
largely  objective  nature  of  the  method  will  reduce  the  inconsistency  between  forecasters 
and  lessen  the  impact  of  inexperience  on  the  reliability  of  the  product. 

Although  this  method  was  developed  for  turbulence,  similar  techniques  using  a 
composite  chart  and  layer  viewer  may  provide  additional  benefits  in  other  areas  of 
forecasting  such  as  for  aircraft  icing.  Additionally,  the  layer  viewer  can  be  used  to  aid 
the  forecaster  in  synoptic  scale  analysis  and  forecasting. 

5.2  Recommendations 

5.2.1  Recommendations  to  AFWA  and  15^^  OFFS.  This  method  appears  to  offer  a 
superior,  comprehensive  tool  in  developing  upper  level  turbulence  charts.  It  allows  a 
large  portion  of  the  globe  to  be  quickly  analyzed  for  turbulence.  Furthermore,  a  method 
very  similar  to  the  test  method  has  been  used  operationally  at  the  17*  OWS  for  over  two 
years.  The  OPS  II  script  employed  by  the  17*  OWS  generates  a  composite  chart  and 
layer  viewer  charts  very  similar  to  those  developed  for  this  study.  However,  the  script 
used  at  the  17*  OWS  does  not  include  mountain  wave  turbulence  and  only  used  output 
from  two  adjacent  layers  instead  of  the  three  layers  used  in  this  study. 
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It  is  strongly  recommended  AFWA  run  the  mountain  wave  forecast  model  VI.  1 
and  V2.1  for  the  troposphere  and  converts  that  output  into  a  format  suitable  for  use  by 
OPS  II.  If  it  is  necessary  to  make  the  data  match  gridded  output,  it  is  suggested  AFWA 
use  the  same  approach  employed  in  this  study.  Furthermore,  the  layer  viewer  html  and 
javascript  code  should  be  provided  to  the  operational  weather  squadrons  for  their 
immediate  implementation.  The  code  developed  for  this  method  can  be  easily  modified 
by  the  weather  squadron  web  masters  and  used  to  supplement  or  replace  current 
techniques. 

Finally,  to  better  meet  the  needs  of  aircrews,  it  is  recommended  that  weather 
packages  include  the  traditional  turbulence  chart,  as  generated  by  the  test  method,  along 
with  several  layer-specific  charts  based  on  unanalyzed  model  output  near  the  flight  level 
of  the  aircraft.  This  provides  the  pilot  with  a  general  idea  of  the  turbulent  regions,  as  well 
as,  the  areas  of  potential  turbulence  at  the  levels  they  are  flying.  Providing  both  charts 
will  offer  the  most  complete  picture  of  where  the  aircrew  might  experience  turbulence 
and  what  their  best  course  of  action  might  be  to  avoid  the  turbulence  if  necessary. 

5.2.2  Future  Research  Recommendations.  There  are  numerous  additional  research 
studies  that  may  be  performed  to  further  test  the  methods  developed  here.  The  test 
performed  in  Table  7  suggests  the  indices’  thresholds  may  vary  based  on  levels  of  the 
atmosphere  and  geographical  locations  of  the  aircraft.  Additional  studies  to  determine  if 
the  thresholds  vary  based  on  those  criteria  might  be  performed.  Furthermore,  this  study 
used  GFS  model  data  to  determine  the  Ellrod-2  index  and  initialize  the  mountain  wave 
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forecast  model.  Using  higher  resolution  model  data  might  provide  more  accurate  results 
than  using  the  data  chosen  here. 

This  study  focused  on  determining  regions  of  moderate  turbulence  using  data  over 
a  five  month  period  to  calibrate  the  indices.  Additionally,  the  study  only  verified 
forecasts  for  two  weeks  in  the  fall.  This  limited  study  was  adequate  to  demonstrate  the 
relative  skill  of  the  test  method  to  operational  forecasts  but  further  studies  using  a  larger 
dataset  might  allow  accurate  calibration  of  severe  turbulence  thresholds  and  it  might 
provide  insight  into  how  the  index  thresholds  change  with  different  seasons.  Future 
studies  might  also  use  PIREPS  that  include  aircraft  type  to  more  accurately  assess 
intensity  reports.  Also,  no  steps  were  taken  to  initialize  or  verify  the  GFS  model  using  to 
compute  the  indices  before  making  the  forecasts.  Taking  those  steps  might  further 
increase  the  accuracy  of  the  forecast  method. 

This  study  only  tested  three  indices  in  predicting  CAT,  using  the  Ellrod-2 
turbulence  index  as  a  means  of  determining  CAT  caused  by  Kelvin-Helmholtz 
instabilities.  There  are  many  other  indices  developed  to  predict  CAT  caused  by  Kelvin- 
Helmholtz  instabilities.  Euture  studies  might  consider  using  a  different  index  or  using  an 
ensemble  forecast  of  sorts  with  several  different  indices  to  determine  if  they  provide 
more  accurate  results.  These  ensembles  might  rely  on  several  indices  on  a  layer 
indicating  turbulence  as  opposed  to  this  study’s  reliance  on  several  different  layers 
depicting  turbulence  to  mark  an  area  with  turbulence.  Additionally,  future  versions  of  the 
MWEM  or  other  mountain  wave  models  might  become  available  that  would  require 
testing  to  determine  their  use. 
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Finally,  this  method  is  a  step  towards  completely  automating  the  CAT  product 
produced  by  weather  squadrons.  It  seems  a  similar  method  using  sets  of  rules  could  be 
developed  to  produce  a  chart  without  any  forecaster  involvement  at  all.  This  method 
might  reduce  the  forecaster’s  role  to  quality  control  checking  the  chart  and  adjusting 
regions  based  on  model  error  in  initialization  or  verification  purposes  only.  This  would 
further  simplify  the  CAT  forecast  process  and  essentially  eliminate  the  forecaster  from 
being  in  the  loop  for  CAT  prediction  altogether. 
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Appendix  A:  Detailed  Computations 


1.  Computations  involved  in  MWFM  Vl.l  (Bacmeister  et  al.  1993).  This  is  designed  to 
provide  a  basic  understanding  of  the  underlying  equations  involved  in  computing 
turbulence  with  the  Mountain  Wave  Forecast  Model. 


First,  the  local  saturation  limit  is  computed 


^savA^) 


Then,  a  first  guess  of  the  wave  displacement  profile  at  a  level  is  made  using 


p{z)N ,{z)U ^Az) 


nl/2 


<5'*(z  +  Az)  =  Si^{z) 

p{z  +  Az)Ni^  (z  +  Az)Uj_.i^  (z  +  Az) 

If  this  value  is  more  than  the  saturation  limit,  it  is  reduced  to  the  saturation  limit. 

To  compute  the  value  at  the  surface,  use  the  minimum  of  saturation  value  at  the 
surface  and  4a^  where  a^.  is  a  ridge  height  parameter.  Now,  the  momentum  flux  is 
approximated  by 


^,{z)^ap{z)N,{z)U^Az) 


L 


and  this  is  used  to  compute  turbulence  by 


KE 


Turbulence 


Az 


Z  +  - 

V  2  J 

The  definition  of  the  terms: 


"^sat’k 


local  saturation  limit  of  wave-induced  vertical  displacement 


perpendicular  wind  component  to  the  ridge  based  on  height 
Briint-Vaisala  frequency 

Si  first  guess  of  wave-induced  vertical  displacement  at  a  level 

p  density  of  the  air  at  a  level 

a  dimensionless  ridge-shape  factor 

L  horizontal  length  representing  extent  of  wave  disturbance 

average  momentum  flux  at  level  k 

^^Turbuience  ^hc  turbulcncc  produced  by  waves  from  lost  momentum  flux 
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2.  Haversine  method  used  to  ealeulate  distanees  on  this  sphere.  Method  as  presented  in 
Sinnott  (1984)  used  extensively  to  ealeulate  distanees.  Funetion  passes  two  sets  of 
latitude  and  longitudes  and  returns  the  distanee  between  them. 

^lat  ~  ~  ^lon  ~  ^^^2  ~ 


( d,  ^ 

2 

f  d,  ^ 

a  =  sin 

'  '^lat  ' 

1 

+  cos(/a^i )  cos(/a^2 ) 

'  ^lon 

1  2 

1  2  j 

f 

c  =  2  arctan 

V 

dist  =  earth 


Vl-a  j 


radius  x  c 
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Appendix  B:  CONUS  what2do  Script 


This  is  the  what2do  script  used  by  this  research  to  run  the  MWFM  VI .  1  and  V2. 1 . 
The  what2do  files  are  located  in  the  /MWFM/source/jobs/what2do/  directory.  Its  primary 
function  is  to  set  up  the  variables  which  control  how  the  experiment  runs.  To  run  a  new 
experiment,  one  must  modify  several  other  scripts  including; 

1 .  /MWFM/source/jobs/what2do/special_cases.pro:  set  up  a  new  special  case 
with  this  number  by  copying  and  pasting  another  special  case  and 
modifying  it  as  needed. 

2.  /MWFM/source/jobs/progs/at_glider.pro:  add  new  special  case. 

3.  /MWFM/source/jobs/drivers/drv_glider.pro:  add  new  special  case. 

4.  /MWFM/source/core/frcst_mtnwave.pro;  add  new  special  case. 

PRO  WHAT2DO_BELSONTEST,date=date,pole=pole,xrang=xrang  $ 


,field=field,levels=levels,latlonlims=latlonlims  $ 
,what=what,source=source,grid=grid  $ 

,sequence=sequence,special=special  $ 

,psfile=psfile,mwfmversion=mwfmversion  $ 


,srcplane=srcplane,strictl8=strictl8,map_ann=map_ann  $  $ 

,plot_items=plot_items 


;+ 


(  )  (  I 

M  )  W  (  E  )  M  I 

(  )  (  I 

_ n _ A\ _ n _ 

Mountain  Wave  Eorecast  Model 


;  WHAT2DO_BEESONTEST 


;  template  what2do_xxxx.pro  routine  to  set  up  a  new  MWEM  experiment:  replace 
;  "xxxx"  with  your  experiment  name  and  change  various  what.*  values  below  to 
;  tune  to  your  specific  experiment,  then  set  up  your  own  drv_xxxx.pro 
;  (MWEM  1.*)  or  specialcase  in  special  cases.pro  (MWFM  2.*). 
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;  MWFM  USAGE  e.g.,  special_cases.pro  (source/jobs/what2do) 
;  LOCATION  source/jobs/what2do 


;  HISTORY 

;  2002081 1  SDE:  created  as  a  standardized  what2do  template  from  eleaned, 
;  eommented,  and  doeumented  what2do_usarayforeeast.pro 


;  KEYWORDS/PARAMETERS 


;  date 

9 

9 

9 

;  pole 

9 

9 

;  xrang 

9 

9 

9 

9 

!  field 

9 

9 

9 

9 

9 

9 

9 

9 

;  levels 


input  of  standard  eharaeter  date  in  YYMMDD  or  YYMMDDZZ 
format:  gets  plugged  into  what.date.  You  should  always 
pass  a  date,  otherwise  you'll  get  an  unsatsifaetory 
default 

2-element  veetor  of  the  "map  pole"  or  eenter  of  the  map, 
in  the  format  [longitude,latitude].  If  not  input,  uses  a 
default  value  set  below:  gets  plugged  into  what.pole 
2-element  x  range  veetor  for  GSEC  map  sizing,  used 
explieitly  only  for  map_souroe=0.  Eor  new  IDE  maps  it 
is  only  used  unless  the  (mueh  preferred)  explieit  sizing 
option  "platlonlims"  in  the  pltree  structure  (i.e., 
pltree.platlonlims)  is  not  set.  If  not  input,  uses  a 
default  value  set  below:  gets  plugged  into  what.xrang 
multi-element  veetor  (#  elements  eorresponding  to  the 
number  of  requested  foreeast  plots),  eontaining  integer 
MWEM  wave  field  indieies  eorresponding  to  various  MWFM 
output  fields  that  ean  be  ealeulated.  Type 
FIELD  STRINGS,  /LIST  within  MWEM  to  see  eurrently 
supported  list  of  permitted  wave  field  indiees  and  what 
they  eorrespond  to.  If  not  input,  field  is  set  to 
default  values  set  below:  field  gets  plugged  into 
what,  field 

a  multi-element  veetor  of  exaetly  the  same  size  as  field 
that  represents  the  integer  index  values  of  the  pressure 
at  whieh  you  wish  to  do  this  foreeast.  This  is  dependent 
on  the  pressure  level  array  of  your  partieular  atmospherie 
data  souree:  MWEM  has  pressures  that  go  from  largest 
(low  altitude)  to  smallest  (highest  altitude),  so  that  if 
the  pressure  level  array  is  [1000.0,925.0,850.0,700.0....] 
then  level=[0,l,3]  tells  MWEM  to  plot  3  foreeasts:  first 
on  the  lOOOhPa  surfaee,  seeond  on  the  925hPa  surfaee  and 
the  third  on  the  VOOhPa  surfaee.  The  field  plotted  in 
eaeh  ease  is  given  by  the  field  array  above:  e.g.. 
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;  field=[0,0,0]  will  plot  peak  vertical  displacements  in 

;  each  case.  To  find  out  allowable  pressure  surfaces,  you 

;  will  need  to  read  your  data  in  using  get  glbdata.pro  and 

;  then  print  out  the  returned  plev  array.  If  not  input, 

;  levels  is  set  to  default  values  set  below:  levels  gets 

;  plugged  into  what,  levels 

;  latlonlims  a  4-element  array  that  defines  the  regional 
;  longitude-latitude  bounds  [lonmin,lonmax,latmin,latmax] 

;  within  which  MWFM  will  do  the  forecast.  For  example, 

;  latlonlims=[-10. 0,30. 0,60. 0,70.0]  performs  a  forecast 

;  in  the  longitude  range  lOW  to  3  OF  and  for  latitudes 

;  between  60N  and  VON,  which  corresponds  roughly  to 

;  Scandinavia.  See  get  glbdata.pro  for  more 

;  information/context.  If  not  input,  latlonlims  is  set  to 

;  default  values  set  below:  latlonlims  gets  plugged  into 

;  what,  latlonlims 

;  what  the  output  "what2do"  structure  that  this  routine  returns, 

;  containing  all  the  information  and  swtiches  needed  to 

;  run  this  particular  MWFM  forecast/analysis  run  just  the 

;  way  the  user  wants.  See  the  structure  assignment 

;  statement  below  for  all  the  information  within  it  and 

;  what  it  all  means,  or  else  type  HELP,  what,  /STRUCTURES 

;  after  calling  this  routine  using  what=what 

;  source  specify  the  source  string:  type  GET_GLBDATA,/LIST  to  get 
;  listing  of  currently  supported  atmospheric  data  sources. 

;  e.g.,  source='NMC'  specifies  NCEP/NMC  data.  If  not  input, 

;  source  is  set  to  default  values  set  below:  source  gets 

;  plugged  into  what,  source 

;  grid  the  grid  string  in  GSFC  format:  e.g.,  'GG2%5X2%5'  for  2.5x 
;  2.5  degree  lon-lat  gridding  in  the  file.  Gets  plugged  into 

;  what,  grid 

;  sequence  the  sequence  data  string  in  GSFC  format,  e.g., 'EOF 
;  special  the  special  data  string  in  GSFC  format,  again  identifying 
;  the  specfic  kind  of  atmospheric  data  file  we're  reading 

;  psfile  if  set  true,  activates  postscript  output  via  what.psfile=l 
;  mwfmversion  passes  the  MWFM  version  (1,  2,  ....) 

;  srcplane  string  used  to  signify  the  source  of  aircraft  data  to  be 
;  overplotted  (e.g.,  'ER2') 

;  strict  18  older  keyword  in  which  atmospheric  data  at  18  standard 
;  pressure  levels  is  strictly  enforced.  If  press=-I 

;  then  it  requires  data  at  all  X  levels  (I  think) 

;  map  ann  older  map  annotation  keyword 
;  plot_items  another  plot  annotation  device,  this  time  a  structure 
;  used  to  annote  something  to  a  plot.  For  a  sample,  see 

;  code  in  what2do_solve.pro 
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9 

;AUTHORIZATION  TO  USE  AND  DISTRIBUTE 

9 

;I  hereby  agree  to  the  following  terms  governing  the  use  and  redistribution  of 
;The  NRL  Mountain  Wave  Eorecast  Model  (MWEM),  mountain  wave  forecasting 
;and  display  software,  written  and  developed  by  Stephen  D.  Eckermann  and 
;colleagues  at  Code  7646  at  the  Naval  Research  Eaboratory  in  Washington,  DC. 

9 

;Redistribution  and  use  in  source  and  binary  forms,  with  or  without 
;modification,  are  permitted  provided  that:  (1)  source  code  distributions 
;retain  this  paragraph  in  its  entirety,  (2)  distributions  including  binary 
;code  include  this  paragraph  in  its  entirety  in  the  documentation  or  other 
;materials  provided  with  the  distribution,  and  (3)  all  advertising  materials 
;mentioning  features  or  use  of  this  software  display  the  following 
;acknowledgment:  "This  product  includes  software  written 
;and  developed  by  Stephen  D.  Eckermann  and  colleagues  of  the  Naval 
;Research  Eaboratory  (NRL)."  Neither  the  name  of  NRL  or  its  contributors,  nor 
;any  entity  of  the  United  States  Government  may  be  used  to  endorse  or  promote 
;products  derived  from  this  software,  nor  does  the  inclusion  of  the  NRL  written 
;and  developed  software  directly  or  indirectly  suggest  NRL's  or  the  United 
;States  Government's  endorsement  of  this  product. 

;THIS  SOETWARE  IS  PROVIDED  "AS  IS"  AND  WITHOUT  ANY  EXPRESS  OR 
IMPLIED 

;WARRANTIES,  INCLUDING,  WITHOUT  LIMITATION,  THE  IMPLIED 
WARRANTIES  OE 

;MERCHANT ABILITY  AND  EITNESS  EOR  A  PARTICULAR  PURPOSE. 

9 


PRINT,  '******************************************************' 

wh3.t2cio  bGlsontGst 

PRINT, 

PRINT,  "N*****************************************************' 

IP  (NOT  KEYWORD_SET(mwfmversion))  THEN  mwfmversion=l  ;default  to  old 
MWEM  1.* 

9 

;  we  must  have  a  date 

9 

IP  (NOT  KEYWORD_SET(date))  THEN  BEGIN 
date='00012300' 
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DATE2STRING,  date,  sdate=sdate,  /cen,  hour=hour 

MESSAGE, '  ****=[=**=i==[=  *  WARNING* *********  ',/INEORMATIONAE 

MESSAGE,  'Date  string  not  input  as  keyword  date=date',  /INEORMATIONAL 
MESSAGE,  'Choosing  default  of '+sdate+',  '+hour+'Z',  /INEORMATIONAL 
MESSAGE  '  =i!H==i!H=*H==i=H=*H=g]s^p)  WARNING********** 

',/INEORMATIONAL 

ENDIE 

9 

;  set  the  atmospheric  data  source  as  defined  in  the  GSEC-like  method  by  4 
;  character  strings:  source=major  source  of  the  data,  sequence=basic  type 
;  of  forecast/analysis  from  this  source,  special  is  a  special  identifier 
;  associated  with  sequence  variations,  and  grid  is  the  latitude-longitude 
;  resolution  of  the  atmospheric  gridding  (e.g.,  grid='GG2%5X2'  means  2.5 
;  degrees  longitude  by  2  degrees  latitude) 

9 

IE  (NOT  KEYWORD_SET(source))  THEN  source='NMC'  ELSE  $ 
source=STRCOMPRESS(STRUPCASE(source),/REMOVE_ALL) 

IP  (NOT  KEYWORD_SET(grid))  THEN  grid='GGlXr 
IP  (NOT  KEYWORD_SET(sequence))  THEN  sequence='pre' 

IP  (NOT  KEYWORD_SET(special))  THEN  special=" 


;  check/set  up  field  and  level  arrays  for  all  the  necessary  plots 
;  see  PIELD  STRINGS,  /LIST  for  permitted  fields  and  their  field  code 
;  see  plevs  from  get  glbdata.pro  for  range  of  permitted  level  indices 


IP  NOT  KEYWORD_SET(field)  THEN  BEGIN 
CASE  source  OP 
'ASM':  field=[0,2,0,2,0,2] 

ELSE:  field=[2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2]  ;  Output  16  turbulence  charts 

ENDCASE 
ENDIE 

;  set  pressure  levels  to  output 
;  10=70mb 

;  ll=50mb  12=30mb  13=10mb  14=5mb  15=2mb  16=lmb,  17=0. 4mb 


IP  NOT  KEYWORD_SET(levels)  THEN  BEGIN 
CASE  source  OP 

'ASM':  levels=[  1 0,10,16,16,19,19] 

ELSE:  levels=[5,6,7,8,9,10,l  1,12,13,14,15,16,17,18,19,20]  ;  Output  levels  850-100 
mb 

ENDCASE  ;  array  same  length  as  field 
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ENDIF 


;  check  field  and  level  array  lengths,  then  set  number  of  plots  (nplots) 


IF  (N_FFFMFNTS(field)  NF  N_FFFMFNTS(levels))  THEN  BEGIN 
MESSAGE,  ’  =i=h==i=n==i=**n=**WARNING**********  ’,/INFORMATIONAL 

MESSAGE,  'field  &  levels  do  not  have  the  same  #  elements',  $ 

/INFORMATIONAL 
MESSAGE,  $ 

'longer  array  will  be  truncated  to  make  sizes  the  same',  $ 

/INFORMATIONAL 

MESSAGE  '  y^^J^JS^JJs^Q********** 

',/INFORMATIONAL 

nplots=MIN([N_ELEMENTS(field),N_ELEMENTS(levels)]) 

ENDIF  ELSE  $ 
nplots=N_ELEMENTS(field) 

MESSAGE,  STRTRIM(STRING(nplots),2)+'  MWFM  Plots  have  been  requested',  $ 
/INFORMATIONAL 

IF  NOT  KEYWORD_SET(latlonlims)  THEN  BEGIN 
latlonlims=[-I40.,-30.,20.,60.]  ;set  for  CONUS 

ENDIF  ELSE  $  ;  1 00  150 

IF  N_ELEMENTS(latlonlims  NE  4)  THEN  BEGIN 
PRINT,  'latlonlims  =  ',latlonlims 
MESSAGE, '  latlonlims  must  be  a  4-element  vector' 

RETURN 

ENDIF 


;  set  map  sizes  and  center/"pole"  location 


IF  NOT  KEYWORD_SET(xrang)  THEN  BEGIN  ;pltrec.platlonlims  can  also  be  used 
xrang=[-0. 425,0. 425]  ;set  your  own  preferred  default  xrang  here 

xrang=[-0.25,0.25] 

ENDIF 

IF  NOT  KEYWORD_SET(pole)  THEN  BEGIN 
polel=[-130.0,39.] 

;  pole2=[-I00.0,39.]  ;set  your  own  preferred  map  poles  -  the  pole  I 
pole2=polel  ;&  pole2  arrays  allow  the  user  to  choose  alternating 
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;map  zones,  originally  created  for  east  coast  and 
;west  coast  US  zooms  for  USA  MWFM  forecasts 

ENDIF 


;  mirror  these  values  into  pole,  x  range  and  data  range  values  for  each  plot 
;  in  the  range  1 . .  .nplots 


poles  =  FLTARR(2,nplots) 
xrangs  =  FLTARR(2,nplots) 
drangs  =  FLTARR(2,nplots) 


9 

;  nominal  data  ranges  for  various  MWFM  fields:  see  FIELD_STRINGS,  /EIST  for 
;  details  on  the  physical  quantities  and  units 


drangO  =  [  [0.,I000.] 

,[0.,I0.] 

$ 

,  [0.,3.0] 

$ 

,  [0.,20.] 

$ 

,  [0.,20.] 

$ 

,  [0.,2.  ]  ] 

$  ;peak  vertical  displacement  amplitudes  (m) 
;Eliassen  Palm  Eluxes  (Pa) 

;turbulence  intensities  (J/m^3) 

;peak  temperature  amplitude  (K) 

;peak  total  horizontal  vel.  amp.  (m/s) 

;peak  vertical  velocity  amplitude  (m/s) 


;  store  poles  and  ranges  for  each  plot 


EOR  i  =  0,nplots/2-l  DO  BEGIN 
poles(0,i*2)  =polel 
poles(0,(i*2+l))  =  pole2 
ENDEOR 

EOR  I  =  0,NPLOTS-1  DO  BEGIn 
xrangs(0,i)  =  xrang 
drangs(0,i)  =  drangO(*,  field(i)  ) 

ENDEOR 

9 

;  other  plot/data-related  default  setups.... 

IE  NOT  KEYWORD_SET(psfile)  THEN  psfile=I 
IE  NOT  KEYWORD_SET(srcplane)  THEN  srcplane=" 

IE  NOT  KEYWORD_SET(strictI8)  then  strictI8=0 
IE  NOT  KEYWORD_SET(map_ann)  THEN  map_ann=0 
IE  NOT  KEYWORD_SET(plot_items)  THEN  plot_items=0 
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;  create  the  "what"  (what2do)  strueture 


what={ 

date:date 


$ 

$ 


;date  string 

$  ;atmospherie  souree  string 

$  ;grid  string 

$  ;sequence  identifier  for  source 

$  ;speeial  identifier  for  souree 

$  ;use  GRIB  format 

press:[1000.00, 925.000, 850.000, 700.000,  $ 

500.000,400.000,300.000,250.000,200.000,150.000,100.000,70.0000,  $ 
50.0000,30.0000,20.0000,10.0000] ,  $ 


souree:souree, 
grid:  grid, 

sequeneeisequenee, 

speeiahspeeial, 

GRIB:1, 


press:- 1 


$  ;set  to  -1  for  non- 18  levels  (e.g.,  TRMM) 
;eomment  out  for  18  standard  pressures 
;set  to  whatever  for  GRIB  pieking  of  levs 


$ 


$ 


nazimuths:18, 
nkvals:2, 
filter:  0, 
minlz:5.0, 
minlh:  100.0, 
amp_erit:200.0, 
latlonlims :  latlonlims, 
only:4,  $ 

units:",  $ 

label:",  $ 

foreeast:'0',  $ 

zerohourforeeast:  1  , 


nplots:nplots, 
interaet:0, 
autoseale:!  , 
densseahl  , 
eolortable:0  , 
metparam:'z' , 


$  ;#  azimuths  in  ray  foreeasts 

;#  horiz.  wavenumbers  if  rayfest 
;filter  MWFM  data? 

$  ;min.  vert  wavlgth  to  filter  (km) 

$  ;min.  horz  wavglth  to  filter  (km) 

$  ;erit.  vdisp  for  mak  rayini 

$  ;  latlonlims  array 

;orig_wd  tag  for  ridge  database 
;unit  string 
;label  string 
;foreeast  flag 
$  ;for  fh=0  assumes  fest  not  anal 

$  ;number  of  plots  to  do 

$  ;interaotively  get  vals 

$  ;automatie  eolor  bar  seale? 

$  ;  1  ->  dens-sealed  momfiux 

$  ;color  table  to  use 

$  ;overplot  what  form  of  data? 


metlevels:fmdgen(30)*2-i-180. ,  $ 


;use  these  eontour  levels 


met_eharsize:0.2 
snatplot:0  , 
nometdata:0  , 
boxes:  1, 
veetors:0  , 
field:field , 


$  ;eharsize  eontour  met  labels 

$  ;plotNAT  supersaturation 

$  ;no  met  data  on  plot 

$  ;plot  ridges  as  ridge  boxes 

$  ;plot  wind  veetors? 

i  ;field  indiees  to  plot 
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$ 

$ 

$ 


levels:levels ,  $ 

strictl8:strictl8  ,  $ 

batchiO,  $ 

traj:0,  $ 

datadumpil,  $ 

powerav:!,  $ 

xavint:0.5, 
yavint:0.5, 
poles:poles, 

drangs:drangs,  $ 

xrangs:xrangs,  $ 

giffde:0,  $ 

gifname:date+'.gif,  $ 
psfdeipsfde,  $ 

pstag:'temp'+date,  $ 
colorps:!,  $ 

noreverse_pscols :  1 ,  $ 

towebsiteiO,  $ 

webscript:'ToWebSite_Glider' 
hardcopy  :0,  $ 

no_flt_trk:  1 ,  $ 

srcplane:srcplane  } 


;level  indicies  @which  to  plot 
;  strict  18  flag 
;older  batch  mode 
;older  trajectory  dump  option 
;dump  data  offline  or  not 
;power  index  4  offline  average 
;averaging  interval  (longitude) 
;averaging  interval  (latitude) 

;map  poles 

;data  ranges  (for  autosoale=0) 

;x  map  ranges 
;export  a  GIF  file? 

;speeify  gif  fde  name 
;export  a  PS  (postseript)  fde? 
;eharaeter  filename  for  psfile 
;oolor  postseript? 

;forees  nonreversal  of  eolor 
transfer  plots  to  web? 

$  ;web  transfer  seript  to  use 
;print  out  a  hard  eopy 
;do/don't  plot  flight  traek 


RETURN 

END 


70 


Appendix  C:  GrADS  Scripts 


There  were  three  files  neeessary  to  use  GrADS  for  this  researeh.  The  first  defines 
the  metadata  file.  It  is  ealled  either  testl2.etl  or  test24.otl.  The  seeond  uses  the  GrADS 
seripting  language  to  build  the  different  layers  used  in  the  layer  viewer.  The  third  file 
uses  the  GrADS  seripting  language  to  verify  the  foreeasts  by  plotting  PIREPS  on  the 
ehart. 

1. 

DSET  ^outl2.dat 
UNDEE  9999 

TITEE  Knapp  Ellrod  and  Mountain  Wave  Indiees 

XDEE  ill  EINEAR -140.0  1.0 
* 

YDEF41  EINEAR  20.0  I 
* 

ZDEE  14  EINEAR  100  50 

H= 

TDEE  I  EINEAR  06NOV2003  IDY 
* 

VARS  4 

h  14  99  average  height  of  layer 
ke  14  99  Knapp  Ellrod  Index 
mwl  14  99  Mountain  Wave  VI 
mw2  14  99  Mountain  Wave  V2 
ENDVARS 

2. 

*  This  program  builds  the  desired  files  the  way  I  want  it  to 

prompt  'Enter  the  desired  date: ' 

pull  datel 

hr=’00' 

time=12 

date2=datel'_'hr 

'set  annot  O' 

'set  baekground  1' 

'set  xlopts  O' 

'set  ylopts  O' 

'set  line  O' 
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'set  grid  on  5  O' 

'c' 

while  (time  <=  24) 
level=700 
llevel=750 

'open  'date2'/test'time'.etr 

say  '***  Opened  'date2'  test'time'.etl  for  analysis  ***' 
while  (level  >=  100) 

'set  mproj  nps' 

'setmpvals  -120.5  -75.5  25  52' 

'e' 

'set  lev  'level 
'set  map  Oil' 

'set  mpdset  mres' 

'set  gxout  shaded' 

'set  elevs  3  6  8.37' 

'set  eeols  1  7  8  2' 

'd  ke' 

'set  gxout  eontour' 

'set  elab  off 

'set  elevs  0.0014  99.0' 

'set  eeolor  4' 

'd  mwl' 

'set  elab  off 

'set  elevs  0.0752  0.103' 

'set  eeolor  9' 

'd  mw2' 

'draw  title  lnitial_Date='date2',  Fest_time='time'Z,  Layer='lever-'llevermb' 
'draw  xlab  Shaded=Ellrod-2  Blue=MWFM  VI .  1  Purple=MWFM  V2. 1' 

'wi  'date2'/'time'_leveriever.jpg' 
level=level-50 
llevel=llevel-50 
endwhile 
'e' 

'set  lev  750' 

'set  gxout  shaded' 

'set  elevs  1  2' 

'set  elab  off 
'set  eeols  13  2' 

'd  h' 

'draw  title  lnitial_Date='date2',  Fest_time='time'Z,  COMPOSITE  ANAEYSIS' 
'draw  xlab  Yellow=moderate  turbulenee  Red=Severe  turbulenee  (possible)' 
'wi  'date2'/'time'_all.jpg' 

'e' 

'elose  r 
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time=time+12 

endwhile 


3. 

*  This  program  builds  the  desired  files  the  way  I  want  it  to 

*  Variables  to  decipher  PIREP  location 
say  'Hello.  PIREP  display  program.' 
x_len=10 

y_len=7 

x_int=x_len/(  1 40-3 0) 

y_int=y_len/(60-20) 

datel=l 

hr='00' 

*  Sets  up  the  basics  of  the  map 
'set  xsize  880  680' 

'set  mproj  scaled' 

'set  annot  O' 

'set  background  1' 

'set  xlopts  O' 

'set  ylopts  O' 

'set  line  O' 

'set  grid  on  5  O' 

'set  background  1' 

'setrgb  16  225  225  225' 

'setrgb  17  200  200  200' 

'setrgb  18  175  175  175' 

'c' 

level=750 


*  Loops  until  the  quit  option  is  chosen 
while  ((datel  !=  0)  &  (datel  !=  ")) 

say  'PIREP  Verification' 

prompt  'Enter  the  desired  date  (0  signifies  quit): 

pull  datel 
date2=dater_'hr 
if  ((datel  !=  0)  &  (datel  !=")) 
tme=12 

while  (tme  <=  24) 

'set  mproj  scaled' 

'set  annot  O' 

'set  background  1' 

'set  xlopts  O' 

'set  ylopts  O' 
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'set  line  O' 

'set  grid  on  5  O' 

'c' 

dates =date2'/ test'tme' .  etl' 

'open  'dates 

say  '  Opened  'date2'  test'tme'. etl' 

'c' 

'draw  title  Date  'date2'  file=test'tme'.etl' 

'set  lev  'level 
'set  map  15  1  1' 

'set  mpdset  mres' 

'set  clevs  1  2' 

'set  cools  17  18' 

'set  gxout  shaded' 

'd  h' 

rlec=0 

cl=l 

while  (rlec  !=  2) 

date4=date2'/ out'tme' .  txt' 
rl=read(date4) 
rlec=substr(r  1,1,1) 
if  (rlec  !=  2) 
lat=substr(rl,S,5) 
if  (substr(rl,15,l)  =  '  ') 

lon=substr(rl,9,6) 
hour=substr(rl  ,27,2) 
hgt=substr(rl  ,S0,5) 
ti=substr(rl,S6,l) 
else 

lon=substr(rl,9,7) 
hour=substr(rl  ,28,2) 
hgt=substr(rl  ,S  1 ,5) 
ti=substr(rl,S7,l) 
endif 

if  ((lat  >=  25)  &  (lat  <=  50)  &  (Ion  >=-lS0)  &  (lon<=-60)) 
x=(140+lon)*x_int+0.5 
y=(lat-20)  *y_int+0 .75 
if(ti  =  0) 
ti2=0 
ms=l 
endif 
if(ti=l) 
ti2=S 
ms=S 
endif 
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if  (ti  >=  2  &  ti  <=  5) 
ti2=7 
ms=5 
endif 
if  (ti  >=  6) 
ti2=2 
ms=9 
endif 

'set  line  'ti2 

'draw  mark  'ms'  'x'  'y'  0.1' 

'set  string  'ti2 
'set  strsiz  0.2' 

'draw  string  'x+O.Ol'  'y-0.05'  'el 
say  el'  'hour'  'laf  'Ion'  'hgf  'ti 
if(cl=5) 

prompt  'Press  enter  to  see  next  bateh' 

pull  temp 

'e' 

'set  lev  'level 
'set  gxout  eontour' 

'set  mproj  sealed' 

'set  mpdset  mres' 

'set  annot  O' 

'set  baekground  1' 

'set  xlopts  O' 

'set  ylopts  O' 

'set  line  O' 

'set  grid  on  5  O' 

'set  baekground  1' 

'set  map  15  1  1' 

'e' 

'set  elevs  1  2' 

'set  ecolor  3  2' 

'set  elevs  1  2' 

'set  eeols  17  18' 

'set  gxout  shaded' 

'd  h' 

'draw  title  Date  'date2'  file=tesftme'.etr 
el=0 
endif 
el=cl+l 
endif 
endif 
endwhile 

prompt  'Press  enter  to  see  next  time!' 
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pull  temp 
tme=tme+12 
'reinit' 
endwhile 
endif 
endwhile 
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Appendix  D:  Turbulence  Viewer  Code 


The  below  html  and  javascript  code  makes  up  the  layer  viewer.  The  first  segment 
of  this  appendix  lists  the  necessary  directories  and  which  files  go  in  it  and  a  brief 
description  of  the  files  if  the  actual  code  is  not  included.  This  is  followed  by  the  code 
that  makes  up  the  two  main  sections  of  the  framed  script. 


DIRECTORY  STRUCTURE  (d=directory,  otherwise  file) 

Eevels  (d) 

EevellOO.htm  -  basic  html  containing  the  flight  level  associated  with  lOOmb 

i 

Eevel700.htm-  same  as  above  except  associated  with  700mb 
script_output  (d) 

12_levell00.jpg  -  image  file  from  12  hour  forecast  at  100-150  mb  layer 

i 

24_level750.jpg  -  image  file  from  24  hour  forecast  at  100-150  mb  layer 
Scrollbars  (d) 

Scrollbar_0.jpg  -  image  of  scrollbar  with  bar  at  top  to  “animate”  scrollbar 

i 

Scrollbar_6.jpg  -  image  of  scrollbar  with  bar  at  bottom  to  “animate”  scrollbar 
Turb_files  (d) 

Bottom  half  htm  -  sets  up  bottom  frames  for  12  hour  forecast,  code  below 
Bottom_half_24.htm  -  similar  to  Bottom_half.htm  for  24  hour  forecast 
buttonsl.htm  -  engine  that  drives  the  layer  viewer,  code  included  below 
buttonsl_24.htm  -  similar  to  buttonsl.htm  for  24  hour  forecast 
layers.htm  -  sets  up  the  layer  frame  with  initial  images 
layers_24.htm  -  similar  to  layers.htm  for  24  hour  forecast 
levels.htm  -  puts  initial  flight  levels  in  flight  level  frame 
TopTitle.htm  -  puts  a  title  in  the  title  frame 
TopTitle_24.htm  -  similar  to  TopTitle.htm  for  24  hour  forecast 
Turb_Eayer_Viewer.htm  -  main  html  code  that  sets  up  a  top  and  bottom  frame 
Turb_Eayer_Viewer_24.html  -  similar  to  Turb_Eayer_Viewer.htm  for  24  hour  forecast 

Bottomhalfhtm 

<html> 

<frameset  cols="75%,  12%,  1 3%"> 

<frame  src="layers.htm"  name="layers"  scrollmg="yes"> 

<frame  src="levels.htm"  name="levels"  scrollmg="no"> 

<frame  src="buttonsl.htm"  name="buttons"  scrollmg="auto"> 
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</frameset> 

</html> 

buttonsl.htm 

<html> 

<head> 

<title>17  OWS  Multi-layer  evaluator</title> 

</head> 

<script  language=" Javascriptl .  1  "> 

//  (T'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'CC'C'C' 

//  Configurable  serollable  layer  viewer  with  most  manipulations  housed  between  rows  of 
dollar  signs 

//  List  of  the  different  levels  available  and  their  labels 

//To  add  new  levels/layers,  add  the  next  ineremental  pair  and  inerease  LEVELMAX 

level  =  new  Array() 

labels  =  new  Array() 

level[0]="100"; 

labels[0]="530"; 

level[l]="150"; 

labels[l]="450"; 

level[2]="200"; 

labels[2]="360"; 

level[3]="250"; 

labels[3]="320"; 

level[4]="300"; 

labels[4]="280"; 

level[5]="350"; 

labels[5]="250"; 

level[6]="400"; 

labels[6]="220"; 

level[7]="450"; 

labels[7]="200"; 


level[8]="500"; 

labels[8]="180"; 


level[9]="550"; 

labels[9]="160"; 
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level[10]="600"; 

labels[10]="140"; 


level[ll]="650"; 

labels[ll]="120"; 


level[12]="700"; 

labels[12]="100"; 

var  LEVELMAX=12; 

//  Size  of  the  step  to  take  along  scrollbar 
var  SCROEEMAX=6; 
var  IMAGEHEIGHT=494; 

varCHANGEAMT=IMAGEHEIGHT/(SCROEEMAX+2); 
var  ERAMESPERPAGE=4; 

//  Pre-loads  the  array  of  images  for  use  in  scrollbar.  You  can  change  the  folder  path 
//  or  add  steps  for  scrollbar. 

UDScrollImage  =  new  Array() 
for(xy=0;  xy  <=  SCROEEMAX;  xy-H-l-)  { 

UDScrollImage[xy]  =  new  Image() 

UDScrollImage[xy].src  =  "../Scrollbars/Scrollbar_"  +  xy  +  ".jpg" 

} 

//  Pre-loads  the  levels  and  labels.  You  may  have  to  change  the  filename  path 
Turbmaps  =  new  Array() 
for(xy=0;xy  <=  EEVEEMAX;  xy++)  { 

Turbmaps[xy]  =  new  Image() 

Turbmaps[xy].src  =  ".. /scrip  t_output/12_level"  +  level[xy]  +  ".jpg"; 

} 

//  (T'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C'C' 

function  ScrollIt(ElD) 

{ 

//  Variables  that  will  not  be  manipulated 
var  temp,temp2,temp3,temp4; 
var  changelty,changeltx; 
changeltx=0; 
changelty=0; 

temp4=parseInt(document.myPorm.ScrollVal.  value,  10); 
if(UD  ==  1)  { 

changeIty=CHANGEAMT; 
temp4=eval(temp4-l- 1 ); 
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}  else  if  (UD  ==2)  { 

changelty=- 1  *CHANGEAMT ; 
temp4=eval(temp4  - 1 ) ; 

}  else  if  (UD  ==  3)  { 

changeltx=- 1  *CHANGEAMT ; 
temp4=eval(temp4  - 1 ) ; 

}  else  { 

ehangeItx=CHANGEAMT; 
temp4=eval(temp4+ 1 ); 

} 

if  (temp4  <  0)  { 
temp4=0; 

} 

if  (temp4  >  SCROLEMAX)  { 
temp4=SCROEEMAX; 

} 

document. myEorm.  ScrollVal.  value=temp4 ; 
temp=eval(document.myEorm.newY.  value); 
temp2=temp+changelty; 
temp=eval(document.myEorm.newX.value); 
temp3=temp+changeltx; 
if  (temp2  <  0)  { 
temp2=0; 

} 

if  (temp3  <  0)  { 
temp3=0; 

} 

document .  myE  orm.  new  Y .  value=temp2 ; 
document.myEorm.newX.value=temp3; 
for  (var  i=0;  i<  ERAMESPERPAGE;  i++)  { 

parent.layers.frames[i].scroll(temp3,temp2); 

} 

document.UDScroll.src=UDScrollImage[temp4].src; 


function  ScrollJump(xy) 

{ 

//  Jumps  to  that  point  on  the  scrollbar 

var  changepomt=xy*CHANGEAMT; 
for  (var  i=0;  i<ERAMESPERPAGE;  i++)  { 

parent.layers.frames[i].scroll(0,changepomt); 

} 

document.UDScroll.src=UDScrollImage[xy].src; 

document.myEorm.ScrollVal.value=xy; 

document.myEorm.newY.value=changepomt; 
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} 


function  MajorScroll(ud) 

{ 

var  nextone=parseInt(document.myFomi.UDVal.  value); 
if  (ud  ==  1)  { 

nextone— ; 

}  else  { 

nextone++; 

} 

if  (nextone  <  0)  { 

nextone  =  0; 

} 

if  (nextone  >  (LEVELMAX-3))  { 
nextone  =  EEVEEMAX-3; 


document. myE  orm.UD  V  al.  value=nextone; 
var  bullet; 

for  (var  i=0;  i  <  ERAMESPERPAGE;  i++)  { 
bullet=i+nextone; 

with  (parent.layers.frames[i]. document)  { 
open() 

writeln("<img  src="+Turbmaps[bullet].src+">"); 
close(); 

} 

with  (parent.levels.frames[i]. document)  { 
open() 

writeln("<html><font 

size=6><center><B>EE"+labels[bullet]+"</B></center></font></html>"); 

elose(); 

} 

parent.layers.frames[i].scroll(0,parselnt(document.myEorm.newY.  value,  10)); 

} 

} 

</script> 

<form  name="myEorm"  method=post> 

<input  type=hidden  name="newY"  value=0> 

<input  type=hidden  name="newX"  value=0> 

<input  type=hidden  name="ScrollVal"  value=0> 

<input  type=hidden  name="UDVal"  value=0> 

</form> 
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<map  name="NSmap"> 

<area  shape="rect"  coords="0,0,15,14"  onClick="ScrollIt(2)";> 

<area  shape="rect"  coords="0,15,15,33"  onClick="ScrollJump(0)";> 
<area  shape="rect"  coords="0, 34, 15,53"  onClick="ScrollJump(l)";> 
<area  shape="rect"  coords="0,54, 15,73"  onClick="ScrollJump(2)";> 
<area  shape="rect"  coords="0,74, 15,93"  onClick="ScrollJump(3)";> 
<area  shape="rect"  coords="0,94,15,l  13"  onClick="ScrollJump(4)";> 
<area  shape="rect"  coords="0,l  14,15,132"  onClick="ScrollJump(5)";> 
<area  shape="rect"  coords="0,133,15,151"  onClick="ScrollJump(6)";> 
<area  shape="rect"  coords="0,152,15,166"  onClick="Scrolllt(l)";> 
</map> 


<image  name="UD Scroll"  border=l  width="15"  height="166" 
src="../Scrollbars/Scrollbar_0.jpg"  usemap="#NSmap"> 

<center> 

<b>Click  on  scrollbar  to  move  up  and  down  inside  the  layers</b> 
<br><br><br><br><br> 

<form  name="myForm2"  method=post> 

<input  type=button  value="  UP  "  onClick="MajorScroll(l)";> 
<br> 

<input  type=button  value="DOWN"  onClick="MajorScroll(2)";> 
</form> 

<b>Click  on  buttons  to  move  shift  the  viewable  layers  up  or  down</b> 
</html> 
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Appendix  E:  FORTRAN  Code  Used  in  Computing  Indices 


The  below  is  the  FORTRAN  eode  used  to  eompute  the  indiees  and  save  the  data 
in  the  proper  format.  It  eonsists  of  the  funetions  ealled  by  a  main  program  whieh  simply 
reads  in  the  model  data  fdes  and  saves  the  indiees  for  display  using  GrADS. 


MODULE  KEmodule 
CONTAINS 

!  Subroutines  used  with  KE.f90  to  eompute  indiees  version  3 


SUBROUTINE  Haversme(latl  ,lat2,lonI  ,lon2,dist) 

!  Uses  the  Haversine  method  to  eompute  distanee 
IMPLICIT  NONE 

!  Deelarations  passed  from  main  program 
REAE,  INTENT(IN)::  latl,lat2,lonl,lon2 
REAE,  INTENT(OUT);;  dist 

!  Eoeally  used  deelarations 

REAE  a,o,dlon,dlat,ER,PI,latlr,lat2r,lonlr,lon2r 

!  Exeeutable  portion 
ER=637 1000.0 
PI=4*atan(l.) 

!  Converts  degrees  to  radians 
latlr=latl*PI/180 
lat2r=lat2*PI/180 
lonlr=lonl*PI/180 
lon2r=lon2  *  PI/ 1 8  0 

!  Computes  the  differenee  in  latitude/longitude  between  points 

dlat=lat2r-latlr 

dlon=lon2r-lonlr 
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!  Calculates  the  distanee 

a=(sin(dlat/2))**2+cos(latlr)*eos(lat2r)*(sin(dlon/2))**2 

c=2*atan(sqrt(a)/sqrt(  1  -a)) 

dist=ER*c 


END  SUBROUTINE  Haversine 
! _ 


SUBROUTINE  Ellrod2 

(U,V,H, BEAT, ELAT, BLON,EEON, USX, usy, EVES, DY,BXSZ,xsz,KE, HU, KEmax) 

!  This  subroutine  calculates  the  Ellrod-2  index  and  returns  a 
!  final  KE  value.  Unlike  AEWA,  this  index  takes  the  average 
!  of  two  levels  and  converts  that  to  a  "layer"  value 

IMPLICIT  NONE 

!  Declare  calling  parameters 

REAL,  INTENT(IN)::  U(:,:,:),  V(:,;,:),  DY,  BXSZ 

INTEGER,  INTENT(IN)::  BEAT,  ELAT,  BLON,  ETON,  EVES,  usx,  usy,  xsz 

REAL,  INTENT(INOUT)::  KE(:,:,:),  HU(:,:,;),  KEmax(:,:) 

!  Locally  stored  variables 
INTEGER  flag,i,j,k,a,b,lat,lon,dx,x,y,z 

REAL  kemx,dudx,dudy,dvdx,dvdy,dudz,dvdz,dst,dsh,def,div,vws,dist 
REAL,  ALLOCATABLE::  H2(:,:,:) 

!  Allocates  a  temporary  array  for  storing  the  level  heights 

x=SIZE(H,l) 

y=SIZE(H,2) 

z=SIZE(H,3) 

ALLOCATE(H2(x,y,z)) 

!  Calculates  the  Ellrod-2  index  for  every  grid  point 

!  Calculates  for  an  extra  row/eolumn  to  allow  PIREP  interpolation  to  work  on  border 
fiag=0.0 
DO  k=2,LVLS 
a=l 

DOj=BLAT,(ELAT-l),-l 

b=l 

lat=91-j 

DO  i=BLON,(ELON+l) 

IP  (xsz  .EQ.  1 8 1)  THEN  !  To  account  for  two  different 
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!  file  sizes 


lon=i-181 
ELSE 

IE  (i.GT.  181)  THEN 
lon=i-361 
ELSE 
lon=i-l 
END  IE 
END  IE 

CALL  Haversine((l .0*lat),(l .0*lat),(l .O*lon),(l .0*lon+BXSZ),dist) 
dx=2*dist 

dudx=(U((i+ 1  ),j  ,k)-U((i- 1  ),j  ,k))/dx 
dvdx=(V  ((i+ 1  ),j  ,k)- V((i- 1  ),j  ,k))/dx 
dudy=(U(i,G+l),k)-U(i,G-l),k))/DY 
dvdy=(V(i,G+l),k)-V(i,G-l),k))/DY 
dudz=(U(i,j,(k+l))-U(i,j,(k-l)))/(H(i,j,(k+l))-H(i,j,(k-l))) 
dvdz=(VG,j,(k+l))-V(i,j,(k-l)))/(HG,j,(k+l))-HG,j,(k-l))) 
dst=dudx-dvdy 
dsh=dvdx+dudy 
def=SQRT(dsh*  dsh+dst*  dst) 
div=dudx+dvdy 

vws=S  QRT(dudz  *  dudz+dvdz  *  dvdz) 

KE(b,a,(k- 1  ))=vws *(def-div) *  1 0 *  * 7 
IE  (KE(b,a,(k-l))  <  0.0)  THEN 
KE(b,a,(k-l))=0.0 
END  IE 
b=b+l 
END  DO 
a=a+l 
END  DO 
END  DO 

PRINT  *, 'Averaging  height  and  KE  values  to  mateh  MWEM  output  levels' 
DO  k=2,LVLS 
kemx=0.0 
DO  b=l,usx 
DO  a=l,usy 

IE  (KE(b,a,(k-l))  .GT.  kemx)  THEN 
kemx=KE(b,a,(k-l)) 

END  IE 

IE  (KE(b,a,(k-l))  .GT.  KEmax(b,a))  THEN 
KEmax(b,a)=KE(b,a,(k-l)) 

END  IE 

HU(b,a,(k-l))=(H(b,a,k)+H(b,a,(k+l)))/2 
END  DO 
END  DO 
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PRINT  Level  ',(k-l),'  of ',(LVLS-1),'  done... max-, kemx 
END  DO 

END  SUBROUTINE  Ellrod2 

! _ 


SUBROUTINE  Mwfm  grid 

(BEAT,EEAT,BLON,EEON,xsz,mwlat,mwlon,mwflux,raylat,raylon,rayflux,MWSZ,EV 
ES,endvI  ,endv2,MW  I  ,MW2,MWImax,MW2max) 

!  Eays  the  MWEM  data  onto  a  Ixl  degree  grid 
!  Uses  the  maximum  value  in  a  box  as  the  value  of  a  grid  point 
!  To  prevent  smoothing  of  the  data 

IMPEICIT  NONE 

!  Declares  calling  parameters 

INTEGER,  INTENT(IN)::  BLAT,ELAT,BEON,EEON,EVES,xsz,endvI(:),endv2(:) 
REAE,  INTENT(IN);;  mwlat(:,;),mwlon(:,:),mwflux(:,:) 

REAE,  INTENT(IN)::  raylat(:,:),raylon(:,:),rayflux(:,:) 

REAE,  INTENT(IN);:  MWSZ 

REAE,  INTENT(OUT):;  MWI(:,:,:),  MW2(:,;,;),MWImax(:,:),MW2max(:,:) 

!  Eocal  variables 
INTEGER  i,j,k,a,b,lat,lon,q2 

!  Execution  portion 

PRINT  *, 'Converting  MWEM  raw  data  into  gridded  format' 

DO  k=  I, (EVES- 1) 
a=I 

DOj=BLAT,EEAT,-I 

b=l 

lat=9I-j 

DO  i=BEON,EEON 

IE  (xsz  .EQ.  1 8 1)  THEN  !  To  account  for  two  different 

lon=i-I8I  !  file  sizes 

EESE 

IE  (i.GT.  1 8 1)  THEN 
lon=i-36I 
EESE 
lon=i-I 
END  IE 
END  IE 
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IF  ((a  .EQ.  1)  .AND.  (b  .EQ.  1)  .AND.  (k  .EQ.  1))  THEN 
PRINT  *, 'Bottom  left  lat,Eon:  ',lat,',',lon 
END  IE 

!  V 1 . 1  gridding 
MWl(b,a,k)=0.0 
DO  q2=l,endvl(k) 

IE  ((mwlat(q2,k)  .GE.  (lat-MWSZ))  .AND.  (mwlat(q2,k)  .EE.  (lat+MWSZ))) 

THEN 

IE  ((mwlon(q2,k)  .GE.  (lon-MWSZ))  .AND.  (mwlon(q2,k)  .EE. 
(lon+MWSZ)))  THEN 

IE  (MWl(b,a,k)  .ET.  mwflux(q2,k))  THEN 
MW  1  (b,a,k)=mwflux(q2,k) 

END  IE 
END  IE 
END  IE 
END  DO 

IE  (MWl(b,a,k)  .GT.  MWlmax(b,a))  THEN 
MW  1  max(b,a)=M  W 1  (b  ,a,k) 

END  IE 

!  V2.1  gridding 
MW2(b,a,k)=0.0 
DO  q2=l,endv2(k) 

IE  ((raylat(q2,k)  .GE.  (lat-MWSZ))  .AND.  (raylat(q2,k)  .EE.  (lat+MWSZ))) 

THEN 

IE  ((raylon(q2,k)  .GE.  (lon-MWSZ))  .AND.  (raylon(q2,k)  .EE. 
(lon+MWSZ)))  THEN 

IE  (MW2(b,a,k)  .ET.  rayflux(q2,k))  THEN 
MW2(b,a,k)=rayflux(q2,k) 

END  IE 
END  IE 
END  IE 
END  DO 

IE  (MW2(b,a,k)  .GT.  MW2max(b,a))  THEN 
MW2max(b,a)=MW2(b,a,k) 

END  IE 
b=b+l 
END  DO 
a=a+l 
END  DO 

PRINT  'Level  ',k,'  of , (EVES- 1),'  Done!' 

END  DO 

END  SUBROUTINE  Mwfm  grid 
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SUBROUTINE 

PR_int(MPR, BEAT, BLON, EVES, METHOD, IX, endvl,endv2,BXSZ2,KEmax,MWlmax 
,MW2max,KE,HU,MWI  ,MW2,H) 

!  Interpolates  INDICES  to  PIREPS 
!  Uses  nearest  neighbor  for  Ellrod-2  index 
!  Uses  one  of  two  methods  for  MWFM 
!  METHOD=  1 .  Uses  nearest  neighbor 

!  METHOD=2.  Uses  elosest  point  within  a  small  box 

IMPEICIT  NONE 

!  Deelarations  passed  from  main  program 

INTEGER,  INTENT(IN)::  MPR,BEAT,BLON,LVLS,METHOD,IX,endvl(:),endv2(:) 
REAL,  INTENT(IN):;  BXSZ2,KEmax(:,:),MWlmax(:,;),MW2max(:,:) 

REAL,  INTENT(IN);:KE(:,:,:),HU(:,:,:),MW1(:,:,:),MW2(:,;,:),H(:,:,:) 

!  Deelarations  used  by  this  program 
INTEGER 

i,numpir,Pyear(MPR),Pmon(MPR),Pday(MPR),Phr(MPR),Palt(MPR),Pturb(MPR) 
INTEGER  elvlo,elvlu,mlvlu,mlvlo,gridx,gridy,ni,nj,k,xsz,ysz 
REAL  eoff,Plat(MPR),Plon(MPR),Plvl,lon2,lat2,dist 
REAL 

Pke_l,Pmwl_l,Pmw2_l,dl  ,d2,d3,d4,ewl  ,ew2,mwwl  ,mww2,wtot,tl  ,t2,mld,m2d,Pke_m, 
Pmw  l_m,Pmw2_m 
CHARACTER  (LEN=2)  fcst  hr 
CHARACTER  (LEN=10)  outfn 

!  Executable 

xsz=SIZE(MWl,l) 

ysz=SIZE(MWl,2) 

PRINT  *, 'Reading  the  PIREP  data  file' 

1 0 1  FORMAT(F5 .2,X,F7.2,X,I4,X,I2,X,I2,X,I2,X,I5  ,X,1 1 ) 

OPEN  (UNIT=20,  FILE='../../PIREPS/PIREP.txf,  STATUS='old') 
i=l 

eoff=2.3 

990  IF  (eoff.EQ.  11.11)  GOTO  991 

READ  (20,101)  Plat(i),  Plon(i),  Pyear(i),  Pmon(i),  Pday(i),  Phr(i),  Palt(i),  Pturb(i) 

eoff=Plat(i) 

i=i+l 

GO  TO  990 
991  CLOSE(20) 
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numpir=i-2 

PRINT  *, 'There  are  ',numpir,'  PIREPS  available  in  the  3  hour  window' 

PRINT  *, 'Storing  the  indice  values  for  each  PIREP  over  the  ',(EVES-1),'  levels.' 

PRINT  Approximate  Eevel  height  range  ',HU(I,I,I),'  to  ',HU(I,I,(LVES-I)),'  meters.' 
OPEN  (UNIT=1I,  FIEE='outfde.tmp',  STATUS='old') 

READ(I  I,*)  fcst  hr 
CEOSE(II) 
outfn='out'//fcst_hr//'  .txt' 

OPEN  (UNIT=20,  FILE=outfn,  STATUS='replace') 

102 

FORMAT(A2,X,F6.2,X,F7.2,X,I2,I2,2X,I2,3X,I5,3X,II,4X,F6.2,X,F8.5,X,F8.5,X,F8.5, 

X,F8.I,X,F8.5,X,F8.I) 

103 

FORMAT(A2,X,F6.2,X,F7.2,X,I2,I2,2X,I2,3X,I5,3X,II,4X,F6.2,X,F8.5,X,F8.5,2X,F6.2, 

X,F8.5,X,F8.5) 

104 

FORMAT(A2,X,F6.2,X,F7.2,X,I2,I2,2X,I2,3X,I5,3X,II,4X,F6.2,X,F8.5,X,F8.5,X,F8.5, 

X,F8.I) 

105 

FORMAT(A2,X,F6.2,X,F7.2,X,I2,I2,2X,I2,3X,I5,3X,II,4X,F6.2,X,F8.5,X,F8.5,I9X,F8. 

5,X,F8.I) 

PRINT  *, 'PIREP  REPORT' 

PRINT  *, 'Hr  Fat  Eon  Date  Hour  Hght  Turb  KE  MWI  MW2' 

DO  i=l,numpir 

!  Converts  feet  to  meters 
Plvl=I.0*Palt(i)/3 .28084 
lon2=I.O*FEOOR(Plon(i)) 
lat2=I  .0*FEOOR(Plat(i)) 
gridx=nint(lon2) 
gridy=nint(lat2) 

IF  (IX  .EQ.  1 8 1)  THEN 
ni=gridx+I4I 
EESE 

IF  (gridx  .LT.  0)  THEN 
ni=gridx+362-BEON 
EESE 

ni=gridx+2-BEON 
END  IF 
END  IF 

IF  (gridy  .FT.  0)  THEN 
gridy=gridy-I 
END  IF 

nj=gridy-90+BEAT 
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DO  k=2,LVLS 

IF  ((H(ni,nj,k)  .GT.  Plvl)  .OR.  (k  .EQ.  EVES))  THEN 
IE  (k  .eq.  EVES)  THEN 
elvlu=EVES-I 
elvlo=EVES-I 
mlvlu=elvlu 
mlvlo=elvlo 
ewl=1.0 
mwwl=1.0 
GOTO  805 

EESE  If(K  .EQ.  2)  THEN 
elvlu=I 
elvlo=I 
mlvlu=elvlu 
mlvlo=elvlo 
ewl=1.0 
mwwl=1.0 
GOTO  805 
EESE 
elvlu=k-I 
elvlo=k 
mlvlu=elvlu 
mlvlo=elvlo 

ew  I = I  -(Plvl-H(ni,nj  ,elvlu))/ (H(ni,nj  ,elvlo)-H(ni,nj  ,elvlu)) 
mww  I = I  -(Plvl-HU(ni,nj  ,mlvlu))/ (HU(ni,nj  ,mlvlo)-HU(ni,nj  ,mlvlu)) 
GOTO  805 
END  IE 

EESE  IE  (H(ni,nj,k)  .EQ.  Plvl)  THEN 
elvlu=k-l 
elvlo=k-l 
mlvlu=elvlu 
mlvlo=elvlo 
ewl=1.0 
mwwl=1.0 

PRINT  *,'PIREP  on  one  level' 

GOTO  805 
END  IE 
END  DO 

805  IE  (ewl  .ET.  0.5)  THEN 
elvlu=elvlu 
EESE 

elvlu=elvlo 
END  IE 

IE  (mwwl  .ET.  0.5)  THEN 
mlvlu=mlvlu 
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ELSE 

mlvlu=mlvlo 
END  IE 

CALL  Haversine(Plat(i)  ,lat2  ,Plon(i)  ,lon2 ,  dist) 
dl=dist 

CALL  Haversine(Plat(i),lat2,Plon(i),(lon2+l),dist) 
d2=dist 

CALL  Haversine(Plat(i),(lat2+l),Plon(i),lon2,dist) 
d3=dist 

CALL  Haversine(Plat(i),(lat2+l),Plon(i),(lon2+l),dist) 
d4=dist 

IL  (dl  .EQ.  MIN(dl,d2,d3,d4))  THEN 
ni=ni 
nj=nj 

ELSE  IL  (d2  .EQ.  MIN(dI,d2,d3,d4))  THEN 
ni=ni+I 

ELSE  IL  (d3  .EQ.  MIN(dI,d2,d3,d4))  THEN 
ni=ni 
nj=nj+I 
ELSE 
ni=ni+I 
ni=ni+I 
END  IE 

Pke_l=KE(ni,nj  ,elvlu) 

Pke_m=KEmax(ni,nj) 

Pmw  I_1=MW  I  (ni,nj  ,mlvlu) 

Pmw  I_m=MW  I  max(ni,nj) 

Pmw2_l=MW2(ni,nj  ,mlvlu) 

Pmw2_m=M  W2max(ni,nj  ) 

WRITE  (20,103) 

fcst_hr,Plat(i),Plon(i),Pmon(i),Pday(i),Phr(i),Palt(i),Pturb(i),Pke_l,PmwI_l,Pmw2_l,Pke 
m,Pmw  I_m,Pmw2_m 
WRITE  (*,103) 

fcst_hr,Plat(i),Plon(i),Pmon(i),Pday(i),Phr(i),Palt(i),Pturb(i),Pke_l,PmwI_l,Pmw2_l,Pke 
m,Pmw  I_m,Pmw2_m 
END  DO 
CLOSE  (20) 

END  SUBROUTINE  PR  int 

! _ 
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SUBROUTINE 

Comp_make(COMPMETHOD, USX, usy, EVES, KEL,KEM,KES,MW1L,MW1M,MW1S, 
MW2L,MW2M,MW2S,KE,MWl,MW2,KEmax,MWlmax,MW2max,HF) 

!  Makes  the  composite  index  grid  point  values 
IMPEICIT  NONE 


!  Declarations  passed  from  main  program 

INTEGER,  INTENT(IN)::  COMPMETHOD,usx,usy,LVLS 

REAL,  INTENT(IN)::  KEL,KEM,KES,MW1L,MW1M,MW1S,MW2L,MW2M,MW2S 
REAL,  INTENT(IN):: 

KE(:,:,;),MWl(:,:,:),MW2(:,:,:),KEmax(:,:),MWlmax(:,:),MW2max(:,:) 

REAL,  INTENT(OUT);:  HF(:,:) 


!  Declarations  used  by  this  program 
INTEGER  i,j,k 

PRINT  *,'Comp_make',KEL,KEM,KES,MWlL,MWlM,MWlS,MW2L,MW2M,MW2S 
HF(1  :usx,  1  ;usy)=0.0 
IF  (COMPMETHOD  .EQ.  1)  THEN 
PRINT  *, 'Method  V 
DO  k=2,(LVLS-2) 

DO  i=l,usx 
DO  j=l,usy 

IF  (HF(i,j)  XT.  1.0)  THEN 

IF  ((KE(i,j,(k-l))  .GE.  KEL)  .AND.  (KE(i,j,k)  .GE.  KEL)  .AND. 
(KE(i,j,(k+l))  .GE.  KEL))  THEN 
HF(i,j)=1.0 
END  IF 

IF  ((MWl(i,j,(k-l))  .GE.  MWIL)  .AND.  (MWl(i,j,k)  .GE.  MWIL)  .AND. 
(MWl(i,j,(k+l))  .GE.  MWIL))  THEN 
HF(i,j)=1.0 
END  IF 

IF  ((MW2(i,j,(k-l))  .GE.  MW2L)  .AND.  (MW2(i,j,k)  .GE.  MW2L)  .AND. 
(MW2(i,j,(k+l))  .GE.  MW2L))  THEN 
HF(i,j)=1.0 
END  IF 
END  IF 

IF  (HF(i,j)  .LT.  2.0)  THEN 

IF  ((KE(i,j,(k-l))  .GE.  KEM)  .AND.  (KE(i,j,k)  .GE.  KEM)  .AND. 
(KE(i,j,(k+l))  .GE.  KEM))  THEN 
HF(i,j)=2.0 
END  IF 

IF  ((MWl(i,j,(k-l))  .GE.  MWIM)  .AND.  (MWl(i,j,k)  .GE.  MWIM)  .AND. 
(MWl(i,j,(k+l))  .GE.  MWIM))  THEN 
HF(i,j)=2.0 
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END  IF 

IF  ((MW2(i,j,(k-I))  .GE.  MW2M)  .AND.  (MW2(i,j,k)  .GE.  MW2M)  .AND. 
(MW2(i,j,(k+I))  .GE.  MW2M))  THEN 
HF(i,j)=2.0 
END  IF 
END  IF 

IF  (HF(i,j)  .FT.  3.0)  THEN 

IF  ((KE(i,j,(k-I))  .GE.  KES)  .AND.  (KE(i,j,k)  .GE.  KES)  .AND. 
(KE(i,j,(k+I))  .GE.  KES))  THEN 
HF(i,j)=3.0 
END  IF 

IF  ((MWI(i,j,(k-I))  .GE.  MWIS)  .AND.  (MWI(i,j,k)  .GE.  MWIS)  .AND. 
(MWI(i,j,(k+I))  .GE.  MWIS))  THEN 
HF(i,j)=3.0 
END  IF 

IF  ((MW2(i,j,(k-I))  .GE.  MW2S)  .AND.  (MW2(i,j,k)  .GE.  MW2S)  .AND. 
(MW2(i,j,(k+I))  .GE.  MW2S))  THEN 
HF(i,j)=3.0 
END  IF 
END  IF 
END  DO 
END  DO 
END  DO 
EESE 

PRINT  *, 'Method  2' 


Do  i=I,usx 
DO  j=I,usy 
HF(i,j)=0.0 

IF  (KEmax(i,j)  .GE.  KEE)  THEN 
HF(i,j)=I.O 
END  IF 


IF  (KEmax(i,j)  .GE.  KEM)  THEN 
HF(i,j)=2.0 
END  IF 

IF  (KEmax(i,j)  .GE.  KES)  THEN 
HF(i,j)=3.0 
END  IF 
END  DO 
END  DO 
END  IF 

END  SUBROUTINE  Comp_make 


END  MODUEE  KEmodule 
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Glossary 


AFCCC 

Air  Boce  Combat  Climatology  Center 

AFGWC 

Air  Boree  Global  Weather  Central  (eurrently  AFWA) 

AFWA 

Air  Foree  Weather  Ageney  (formerly  AFGWC) 

AVAR 

Automated  vertical  accelerometer 

AVN 

Aviation  model 

CAT 

Clear-air  Turbulence 

CONUS 

Continental  United  States 

CVG 

Convergence 

CSI 

Critical  Skill  Index 

DBF 

Deformation  shear  to  include  shearing  and  stretching  deformation 

DSH 

Shearing  deformation 

DST 

Stretching  deformation 

FAA 

Federal  Aviation  Administration 

FAR 

False  alarm  rate 

FL 

Flight  Bevel 

FTP 

file  transfer  protocol 

FORTRAN 

Bormula  Translation  programming  language 

GFS 

Global  Borecast  System 

GrADS 

Gridded  Analysis  and  Display  System 

GRIB 

Gridded  Binary  format 

GSM 

Global  spectral  model 

GSS 

Gilbert  Skill  Score 

HiCAT 

High  altitude  clear-air  turbulence  software 

HIRAS 

High-Resolution  Analysis  System 

HSS 

Heidke  Skill  Score 

IDL 

Interactive  Data  Banguage 

ITFA 

Integrated  turbulence  forecast  algorithm 

JPEG 

Joint  Photographic  Experts  Group 

MB 

MegaBytes 

MOG 

Moderate  or  greater 

MSB 

Mean  sea  level 

MWBM 

Mountain  wave  forecast  model 

NCAR 

National  Center  for  Atmospheric  Research 

NCEP 

National  Center  for  Environmental  Prediction 

NGM 

Nested  grid  model 

NMC 

National  Meteorological  Center 

NRB 

Naval  Research  Baboratory 

OPS  II 

Operational  Production  System  version  II 

ows 

Operational  Weather  Squadron 

PIREP 

Pilot  report 

POD 

Probability  of  detection 

POD„ 

Probability  of  detection  of  regions  without  turbulence 

PODy 

Probability  of  detection  of  regions  with  turbulence 

RAOB 

Rawinsonde  observation  program 

Ri 

Richardson  number 

Th 

Ellrod  index  1 

TI, 

Ellrod  index  2 

TSS 

True  skill  score 

USAF 

Elnited  States  Air  Eorce 

VE 

Volume  efficiency 

VWS 

Vertical  wind  shear 

WMO 

World  Meteorological  Organization 
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